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1.  Objective 

a)  Overall  objective  is  to  acquire  fundamental 
understanding  of  a  phenomenon  characterized  by  violent 
fluctuation  of  swirling  flow,  which  is  often  found  to 
occur  in  various  aircraft  engine  components .  This  flow 
instability,  dubbed  here  as  "Vortex  Whistle",  is  one  of 
the  most  subtle  and  treacherous  flow-induced  vibration 
problems  in  gas  turbines.  In  contrast  to  the  other 
well-known  unsteady  flow  problems  in  turbomachinery  such 
as  rotating  stall,  surge,  aeroacoustic  noise  and  flutter, 
at  present,  little  is  known  about  this  phenomenon  -- 
despite  its  importance  to  the  aircraft  engine  overall 
structural  integrity.  The  resultant  vibration  induced 
by  the  "Vortex  Whistle"  can  sometimes  become  so  violent 
that  the  bladings  and  the  structural  members  of  gas  tur¬ 
bines  suffer  serious  damage.  Perhaps  for  the  reason 
that  the  phenomena  have  appeared  in  seemingly  unrelated 
incidents  concealed  under  various  disguises,  so  far  no 
investigations  appear  to  have  been  carried  upon. 

In  the  present  effort,  we  will  conduct  a  comprehen¬ 
sive  and  systematic  investigation  into  the  "Vortex 
Whistle"  with  its  objective  to  offer  a  unifying  explana¬ 
tion  for  this  least  understood  problem  and  to  contribute 
to  assuring  adequate  design  margins  in  order  to  alleviate 
this  severe  flow-induced  vibration  problem  encountered 
in  gas  turbines .  The  entire  program  is  comprised  of 
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theoretical  and  experimental  investigations.  In  Phase  I, 
covering  the  two-year  period  between  April  1,  1978,  and 
March,  1980,  a  theoretical  investigation  has  been  conducted 
and  completed.  Based  upon  this  framework,  we  shall  carry 
out  the  experimental  program  in  the  second  phase  starting 
April  1,  1980.  In  order  to  accelerate  the  pace  of  the 
entire  investigation,  a  part  of  the  experimental  investi¬ 
gation  corresponding  to  the  second  phase  has  been  conducted 
concurrently  with  Phase  I. 

b)  The  objectives  of  the  first  year;  Phase  I  -  1, 

(April  1,  1978  to  March  31,  1979),  the  activity  of  which 
was  reported  in  the  first  annual  report,  were  to  lay  out 
the  basic  analytical  formulation  and  from  it  to  obtain  the 
preliminary  theoretical  explanation  of  the  problem. 

c)  The  objectives  of  the  second  year;  Phase  1-2, 

(April  1,  1979  to  March  31,  1980),  whose  outcome  is  sum¬ 
marized  herein,  are  to  refine  the  foregoing  analysis  and 
resolve  the  several  important  issues  raised  in  the  first 
year.  Throughout  the  first  and  the  second  year,  the  theoreti 
cal  effort  has  been  focused  on  devising  a  flow  model  which 

is  simple  enough  to  be  amenable  to  analysis,  but  still  cap¬ 
tures  the  essential  physics  and  exhibits  the  key  feature 
of  the  "Vortex  Whistle"  phenomenon. 

In  addition,  based  upon  the  results  of  the  analysis, 
the  test  apparatus  has  been  designed,  constructed  and  in¬ 
stalled,  the  data  acquisition  therefrom  comprising  the 
central  part  of  the  second  phase  activity. 


2.  Features  of  "Vortex  Whistle"  Phenomenon 

"Vortex  Whistle"  has  been  known  to  occur  in  various 
gas  turbine  components  such  as  (a)  a  downstream  section 
of  variable  vanes  followed  by  accelerating  flow  (b)  an 
inducer  section  of  centrifugal  compressors  installed  with 
variable  pre-swirl  vanes  and  (c)  turbine  cooling  air 
cavity  where  the  air  enters  through  ports  located  on 
rotating  parts.  The  common  features  of  this  unsteady  flow 
oscillation  may  be  summarized  as  follows: 

(a)  The  most  unmistakable  characteristics  of  "Vortex 
Whistle"  is  that  the  frequency  of  fluctuation  is 
discrete  and  it  becomes  higher  as  the  flow  rate 
increases . 

(b)  It  is  induced  by  high  swirl  flow:  if  the  vane 
angle  is  set  in  such  a  way  as  to  produce  less 
swirl,  whistle  disappears. 

(c)  However,  the  role  of  vanes  in  connection  with 
the  "Vortex  Whistle"  seems  only  to  impart  the 
swirling  motion  to  the  fluid;  in  the  place  of 
vanes,  a  single  or  several  tangential  injection 
of  flow  produces  the  similar  oscillation. 

(d)  The  steady  velocity  distribution  appears  to 
affect  intimately  the  occurrence  of  the  whistle. 
For  example,  small  change  in  duct  configuration 
or  design  change  from  free  vortex  to  forced 
vortex  has  sometimes  succeeded  in  eliminating 
the  pulsation  of  flow. 


(e)  The  amplitude  of  oscillation  often  becomes 


exceedingly  large.  Interestingly  and  curiously 
enough,  in  such  situations  the  steady  flow 
distribution  in  the  radial  direction  --  both 
velocity  and  temperature  --  becomes  markedly 
altered.  For  instance,  for  swirl  within  a 
coannular  passage  between  outer  and  inner  sur¬ 
faces,  the  free  vortex  distribution  of  swirl 
flow  observed  at  points  below  certain  threshold 
swirl  Mach  number  is  found  to  be  converted  into 
a  shape  somewhat  similar  to  a  forced  vortex 
above  the  threshold  Mach  number  --  with  reduced 
swirl  velocity  near  the  inner  surface.  At  the 
same  time,  the  steady  total  temperature  distri¬ 
bution  in  the  radial  direction  exhibits  the 
temperature  drop  as  much  as  30°F  at  the  core  of 
the  vortex;  this  immediately  presents  its  important 
and  intriguing  implications  related  to  Ranque- 
Hilsch  vortex  cooling  effect  (e.g.  Ref.  1). 
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3.  Significant  Achievements  to  Date 

The  following  are  the  significant  results  accrued  in 
the  second  year's  activity  from  April  1,  1979,  to  March  31, 
1980. 

(a)  Analytical  investigation 

The  problem  posed  is  to  study  the  characteristics 
of  disturbance -induced  flow  field  within  an  annular  passage 
between  concentric  circular  cylinders . 

As  reported  in  our  first  annual  report,  in  the 
first  year,  linearized  acoustic  wave  problem  was  analyzed 
and  the  frequency-swirl  relationship  was  obtained;  the  com¬ 
parison  with  the  available  experimental  data  showed  favorable 
quantitative  agreement  with  the  experimental  data  and  repro¬ 
duced  the  trend  stated  in  (a)  of  Section  2.  In  addition,  the 
expression  of  the  steady  streaming  was  derived  using  the 
boundary  layer  approximation  as  a  starting  point;  the  result 
predicted  the  existence  of  the  reversal  in  the  direction  of 
tangential  streaming  and  based  on  this,  we  were  able  to 
explain  the  observed  deformation  of  steady  profile,  referred 
to  Section  2(e). 

Although  these  results  of  the  first  year  were 
highly  encouraging,  the  use  of  the  boundary  layer  approxi¬ 
mation  as  a  starting  point  raised  some  questions  on  the 
expression  of  the  steady  streaming  thus  derived.  The  reason 
is  that  the  streaming  is  obtained  as  a  second  order  quantity 
while  the  boundary  layer  approximation  is  the,  so-called, 
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first  order  approximation  to  the  full  Navier-Stokes 
equations.  Thus,  for  example,  the  effect  of  curvature 
is  neglected  in  the  boundary  layer  approximation  and  we 
have  to  confirm  whether  it  influences  the  steady  stream¬ 
ing  or  not;  the  variation  of  viscosity  due  to  temperature 
fluctuation  has  been  neflected  in  the  first  year's  analysis 
and  this  effect  needs  to  be  assessed. 

To  face  these  problems  once  and  for  all,  in  the 
second  year,  we  decided  to  use  the  Navier-Stokes  equation 
as  the  starting  point  of  the  refined  analysis.  By  using 
the  apparatus  of  a  matched  asymptotic  expansion,  we  derived 
the  expression  of  streaming  anew. 

In  the  main,  this  improved  expression  agrees  es¬ 
sentially  with  the  first  year's  result  based  upon  the 
boundary  layer  approximations,  the  only  major  difference 
being  in  the  presence  of  the  viscosity  fluctuation  in  re¬ 
sponse  to  temperature  unsteadiness.  Even  this  does  not, 
however,  affect  the  values  of  the  threshold  swirl  Mach 
number.  Thus,  we  are  able  to  confirm  the  conclusions  reach¬ 
ed  in  our  initial  efforts,  which  we  summarize  next. 

The  streaming  in  the  tangential  direction  suffers 
a  sudden  reversal  of  its  direction  above  the  threshold  swirl, 
the  physical  reason  being  due  to  the  Doppler  shift  caused 
by  swirl;  and  at  the  same  time,  the  absolute  magnitude  of 
acoustic  streaming  becomes  considerably  increased. 

The  specific  value  of  swirl  at  the  threshold  depends 
on  the  explicit  relationship  between  the  frequency  and  the 
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prescribed  radial  distribution  of  tangential  velocity. 

For  the  steady  free  vortex  distribution  between  co-annular 
cylinders,  which  corresponds  to  the  one  referred  to  (e) 
of  Section  2,  the  threshold  steady  swirl  is  shown  to 
range  from  subsonic  to  supersonic  tangential  Mach  number, 
its  specific  values  being  dependent  upon  the  ratio  of  the 
outer  to  inner  radius  of  the  cylinder  and  the  wave  modes . 

Below  the  threshold,  the  tangential  streaming  on  the  surface 
of  the  inner  cylinder  is  in  the  same  direction  as  the  steady 
swirl;  however,  beyond  this,  the  streaming  abruptly  reverses 
its  direction  and  starts  to  retrogress  in  the  direction 
opposite  to  steady  swirl.  Surely,  then,  this  tends  to 
decrease  the  total  d.c.  component  of  circumferential  velocity, 
which  is  the  sum  of  steady  swirl  with  tangential  streaming, 
its  reduction  being  sizeable  near  the  threshold;  close  to 
the  inner  cylinder,  the  radial  profile  is  converted  into  one 
not  unlike  a  forced  vortex.  This  behavior  appears  to  be 
consistent  with  the  observation  made  in  (e)  of  Section  2. 

We  also  extended  the  analysis  to  include  the  steady 
Rankine  vortex  distribution  within  a  single  pipe,  which 
corresponds  to  the  so-called  Ranque-Hilsch  tube.  There, 
for  the  fundamental  mode  of  disturbance,  any  amount  of  swirl 
always  makes  the  tangential  streaming  near  the  tube  periphery 
rotate  in  the  same  direction  as  the  steady  swirl  itself;  for 
such  a  wave,  no  threshold  swirl  exists.  Its  magnitude  be¬ 
coming  of  considerable  strength,  the  total  d.c.  component 
of  circumferental  velocity  in  a  free  vortex  region  is  in- 
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creased  and  the  entire  Rankine  vortex  is  now  converted 
into  a  forced  vortex. • 

These  conversions  toward  forced  vortex  type  by 
streaming  tend  to  separate  the  flow,  with  the  initially 
uniform  total  temperature,  into  hotter  air  near  the  outer 
radius  and  colder  air  near  the  inner  radius  or  the  center- 
line  --  this  is  the  Ranque-Hilsch  effect. 

The  details  of  the  foregoing  analysis  carried 
out  as  Phase  I  are  written  up  in  a  paper  entitled  "  Steady 
Streaming  in  Swirling  Flow  and  Separation  of  Energy"  to 
be  submitted  for  publication. 


Based  upon  the  effect  of  governing  parameters 
predicted  in  (a),  a  test  rig  has  been  designed  in  order 
to  correlate  the  analysis  with  experiments .  The  primary 
consideration  in  the  design  was  the  establishment  of  a 
well-defined  swirling  flow  in  a  straight  co-annular  tube 
of  sufficient  diameter  to  allow  the  insertion  of  probes 
without  causing  a  major  disturbance  of  the  flow;  the  ratio 
of  the  radii  of  inner  and  outer  tubes  is  to  be  varied  over 
a  range  of  parameters  in  order  to  make  comparisons  with 
the  trend  predicted  by  analysis. 

It  was  decided  that  the  tube  should  be  transpa¬ 
rent  to  allow  flow  visualization  studies  and  a  three- 
inch  outer  diameter  plexiglass  tube  of  30  inches  in  length 
was  selected  (Fig.  1).  The  manifold  section,  into  which 
the  compressed  air  flows ,  is  fabricated  of  transparent 
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plastic.  Eight  nozzles  of  3/8  inches  diameter  equally 
spaced  around  the  circumference  of  the  tube  and  rounded 
at  their  entrance,  direct  this  compressed  air  tangentially 
into  the  main  co-annular  rube.  At  the  other  end  of  the 
30-inch  long  tube,  a  60°  cone-shaped  valve  is  located  to 
regulate  the  amount  of  through  flow. 


Figure  2.  View  of  Vortex  Whistle  Test  Rig  with 
Fixed  Tangential  Inlet. 
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Six  different  inner  tubes  of  varying  radii  are  to  be 
used  in  order  to  investigate  the  effect  of  the  outer/ 
inner  radius  ratio.  The  air  entering  the  manifold  is 
fed  through  a  specially  constructed  acoustic  muffler. 

The  photo  of  the  entire  test  arrangement  is  shown  in 
Fig.  2. 

In  addition  to  the  manifold  of  fixed  inlet 
geometry  described  above,  another  manifold  containing 
variable  guide  vanes  has  been  designed  (Fig.  3);  the 
latter  is  interchangeable  with  the  former,  provides  the 
variation  in  swirl/axial  velocity  ratio,  and  the  bulk 
of  data  is  planned  to  be  taken  with  this  manifold.  The 
variable  guide  vanes  consist  of  24  streamlined  airfoils, 
placed  symmetrically  in  a  circular  array  around  the  side 
plate  of  the  manifold.  The  vane  angles  can  be  set  at  any 
desired  value  between  zero  (purely  radial  flow)  and  65 
degrees,  by  rotating  a  circular  ring  which  contains  24 
small  pivots  fitted  into  slots  provided  in  airfoils .  The 
photo  of  this  variable  vane  manifold  is  shown  in  Figure 
4. 

Currently  check-out  tests  of  ths  test  rig  are 
being  carried  out. 


24  VANES 


Vaneg 


Figure  4.  View  of  Variable  Guide  Vane  Manifold. 
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4.  Implications  of  Conclusions  Obtained  During 

Phase  I  Effort  as  Related  to  Aircraft  Gas  Turbines 

Once  confirmed  by  further  investigations  to  be 
carried  out  in  Phase  II,  the  implications  of  these 
conclusions  obtained  from  Phase  I  program,  as  related  to 
the  aircraft  engine  technology  are  twofold: 

(a)  By  explicit  recognition  of  the  dependence 
of  vortex  whistle  upon  the  governing 
parameters  as  found  in  the  present  in¬ 
vestigation,  it  appears  possible  to 
avoid  the  catastrophic  structural 
failure  by  de-tuning  the  natural  fre¬ 
quency  of  various  engine  components 

from  this  discrete  frequency. 

(b)  The  existence  of  transition  from  free 
vortex  to  non-free  vortex  type  above 
certain  critical  swirl  Mach  number 
implies  that  in  the  steady  aerodynamic 
design  of  rotors/stators,  a  due  considera¬ 
tion  may  have  to  be  given  to  this 
acoustic  streaming. 


5.  Written  Publications 

"Steady  Streaming  in  Swirling  Flow  and  Separation 
of  Energy"  (attached  as  Appendix  and  to  be  submitted 
for  publication) . 


6.  Future  Plans 


(a)  For  remainder  of  present  contrant  (Phase  II, 

April  1,  1980  to  September  30,  1981). 

The  phase  II  program  will  be  comprised  of  (1) 
detailed  data  acquisition  of  'Vortex  Whistle'  by  utilizing 
the  testing  designed  and  constructed  during  Phase  I,  as 
described  in  Section  3(b)  and,  (2)  concurrent  continuation 
and  refinement  of  analysis,  if  necessary. 

The  experimental  program  will  be  initiated  first 
by  installing  the  manifold  of  fixed  geometry  type  and  the 
measurement  of  steady  temperature  and  pressure  distribution 
will  be  carried  out  in  the  test  configuration  where  the 
inner  tube  will  be  removed.  The  compressed  air  with  the 
inlet  pressure  of  10,  15  and  20  psig  will  be  directed  into 
the  manifold  and  flow  field  traverse  will  be  made.  The 
objective  of  this  test  run  is  to  check  out  the  rig  and 
instrumentation  by  making  comparison  with  the  well -documented 
data  taken  in  a  similar  set-up.  Upon  completion  of  this, 
the  manifold  will  be  replaced  by  the  alternative  one  with 
variable  guide  vanes;  this  will  allow  one  to  vary  the  ratio 
of  swirl  velocity  to  axial  velocity.  After  checking  the 
uniformity  of  the  incoming  flow  by  means  of  flow  visualiza¬ 
tion,  we  measure  the  variation  of  steady  velocity  and  tem¬ 
perature  profile  as  the  swirl  Mach  number  is  increased; 
particular  attention  will  be  focused  upon  the  data  acquisi¬ 
tion  near  the  critical  swirl  Mach  number.  In  conjunction 
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with  this,  unsteady  flow  measurement  will  be  carrued  out 
to  define  the  acoustic  characteristics  of  the  vortex 
whistle;  both  the  microphones  placed  externally  and  the 
miniature  kulites  mounted  flush  with  the  internal  surface 
of  the  outer  tube  will  be  used  for  this .  Among  the 
acoustic  signatures  to  be  measured  and  analyzed  are  fre¬ 
quency  of  vortex  whistle,  its  intensity  and  rotating  pat¬ 
terns  of  sound.  In  the  event  where  the  level  of  its  inten¬ 
sity  in  the  test  rig  were  not  sufficiently  high,  it  would 
be  amplified  by  a  speaker,  which  will  provide  an  additional 
source  of  excitement.  Based  upon  the  measured  frequency, 
an  acoustic  suppressor  will  be  designed  and  installed  on 
the  test  rig  and  the  effect  of  sound  suppression  upon  the 
steady  flow  profile  near  the  critical  Mach  number  will  be 
investigated.  As  an  added  related  effort,  the  similar 
effect  upon  the  Ranque-Hilsch  tube  will  be  examined  by 
installment  of  an  acoustic  suppressor  on  the  commercially 
available  vortex  tube,  which  has  been  already  procured  in 
the  Phase  I  period.  With  regard  to  the  analysis,  we  make 
comparisons  of  the  data  with  the  preliminary  analysis 
carried  out  in  Phase  I  and,  if  necessary,  modify  it. 

(b)  Next  increment 

In  the  foregoing  phases;  the  outer  and  inner 
cylinders  simulating  the  outer  and  inner  casing  of  the 
gas  turbine  are  held  stationary.  For  the  actual  aircraft 
engine  turbomachinery ,  the  inner  casing  is  rotating,  of 
course.  This  rotational  effect  is  considered  to  be  im- 


portant  in  the  present  study  of  streaming,  affecting 
the  frequency  parameter.  Thus  in  the  next  increment 
of  the  contract,  the  investigation  of  this  effect  will 
be  proposed. 
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ABSTRACT 


This  paper  concerns  the  steady  streaming  induced  by 
unsteady  disturbances  in  a  swirling  flow  contained 
within  concentric  circular  cylinders  or  a  single  tube. 

The  investigation  is  motivated,  in  the  first  place,  by 
a  newly  observed  phenomenon,  which  reveals  that  the 
acoustic  streaming  is  accountable  for  the  deformation 
of  base,  steady  swirl  profile,  leading  to  the  radial 
separation  of  total  temperature  (the  Ranque-Hilsch 
effect).  This,  in  turn,  offers  a  clue  into  the  hithet- 
to  unheeded  mechanism  of  the  Ranque-Hilsch  effect  itself. 
Starting  from  the  full,  compressible,  unsteady  Navier- 
Stokes  equations,  the  acoustic  streaming  is  studied  by 
the  method  of  matched  asymptotic  expansions;  based  on 
the  results,  experimental  observations  are  explained. 


1.  Introduction 


1 . 1  Background 

The  subject  of  acoustic  streaming  owes  its  origin 
to  Lord  Rayleigh's  landmark  memoir  (1884);  led  by  Faraday's 
observation  (1831)  and  the  patterns  in  the  Kundt  tube,  he  showed 
that  sound  waves  can  and  do  generate  steady  current  through 
the  very  action  of  Reynolds  stresses,  which  are  induced  near 
the  solid  boundary  by  the  periodic  disturbances  themselves. 

For  the  modern  review  of  the  subject  in  general,  we  refer 
to  the  recent  expository  lecture  by  Lighthill  (1978-a) . 

Not  only  can  we  demonstrate  the  acoustic  streaming 
in  assorted  laboratory  experiments  using  a  vibrating  dia¬ 
phragm,  cylinder  and  the  like,  but  it  has  been  suggested  as 
a  possible  explanation  ranging  from  the  roll  torque  effects 
of  rocket  motors  in  flight  (Swithenbank  and  Sotter  1964 ; 

Flandro  1964,  1967)  to  the  blood  flow  phenomena  in  the 
coronary  arteries  (Secomb  1978). 

Of  late,  a  striking  acoustic  streaming  phenomenon, 
with  its  features  alien  to  others,  revealed  itself  unexpect¬ 
edly  in  a  swirling  flow  experiment  (Danforth  1977;  Rakowski, 
Ellis  and  Bankhead  1978;  Rakowski  and  Ellis  1978),  display¬ 
ing  a  grossly  deformed  pattern  of  steady  flow  and  temperature. 
And,  at  the  same  time,  it  afforded  a  glimpse  into  the  dimly 
foreseen  mechanism  of  energy  separation  --  the  Ranque-Hilsch 
effects.  Although  the  phenomenon  was  detected  in  a  test 


2 


rig  called  an  annular  cascade  simulating  flow  in  aircraft 
engines  and,  as  a  matter  of  fact,  cropped  up  in  the  check¬ 
out  tests  as  an  undesirable  side-effect  to  be  eliminated 
later,  its  undeniable  significance  --  in  divulging  the  clues 
to  connect  the  acoustic  streaming  with  thermal  effects  -- 
appears  to  transcend  beyond  the  special  interest  of  turbo¬ 
machinery  technology  and  merit  wider  scientific  attention. 

Here  we  outline  its  layout  briefly.  The  annular  cas¬ 
cade  was  conceived  with  the  objective  to  investigate  some 
aero-elastic  aspects  of  compressor  bladings  in  a  non-rotat¬ 
ing  environment.  As  a  whole,  it  is  in  the  shape  of  a 
stationary,  annular  conduit  formed  between  inner  and  outer 
casings:  first,  air  enters  the  vehicle,  axially  and  uniformly, 
and  is  immediately  imparted  a  tangential  motion  by  passing 
through  variable  swirl  vanes,  the  adjustment  of  whose  angle 
induces  the  change  in  swirl;  then  , it  flows  spiralling  aft, 
through  the  transition  piece,  to  the  test  section  where  an 
array  of  removable  test  airfoils  are  normally  mounted  on  outer 
casings  in  the  circumferential  direction  and  in  a  cascade 
arrangement;  and  finally,  upon  being  realigned  in  the  axial 
direction  by  deswirl  vanes,  the  air  exhausts  to  the  exit. 

During  the  check-out  test  of  the  vehicle,  the  pres¬ 
ence  of  loudly  audible,  unsteady • flow  was  immediately  un¬ 
covered.  The  disturbance  became  manifest  beyond  certain 
conditions  called  an  acoustic  boundary.  It  was  organized, 
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periodic,  and  spinning  circumferentially  with  the  first 
tangential  mode;  the  amplitude  of  total  pressure  exceeded 
20%  of  steady  state  levels,  an  intense  fluctuation  indeed. 

Its  fundamental  frequency  was  in  the  range  of  300  to  400 
Hz,  aerodynamically  ordered  and,  in  fact,  was  found  to 
increase  almost  proportionally  to  the  swirl,  a  point  to 
be  made  here  and  recalled  later.  As  the  swirl  was  increased, 
several  higher  harmonics  were  found  to  accompany  this  funda¬ 
mental  frequency.  The  measurements  were  taken  without  test 
airfoils  installed  in  the  test  section.  Hence, nothing  lay 
in  the  way  of  the  flow  path  between  the  upstream  swirl  vanes 
and  the  downstream  deswirl  vanes .  Before  we  go  further ,  we 
can  not  too  strongly  emphasize  the  fact  that,  despite  its 
aim,  all  the  components  of  the  annular  cascade  are  not  rotat¬ 
ing,  but  stationary. 

Among  the  other  effects  of  this  vigorous  pulsation  in 
swirling  flow,  the  phenomenon  that  arouses  our  attention  is  the 
unexpected  change  of  steady-state  or  time-average  components 
of  the  flow  field, or  its  'H.c.M  parts.  When  the  swirl  was  small 
and  outside  of  the  acoustic  boundary,  the  steady-state  tan¬ 
gential  velocity  distribution  in  the  radial,  direction  was  in 
the  form  of  a  free  vortex,  with  the  obvious  exception  of  thin 
boundary  layers  found  near  the  inner  and  outer  surfaces;  the 
steady-state  tqtal  temperature  was  uniform.  The  former  was 
what  had  precisely  been  intended  in  the  design,  the  latter 
as  expected.  However,  when  the  swirl  was  increased  beyond 


the  acoustic  boundary,  then  above  a  certain  swirl,  the  tan¬ 
gential  velocity  near  the  inner  wall  became  abruptly  re¬ 
duced  to  a  considerable  extent,  the  radial  profile  trans¬ 
figured  from  a  free  vortex  into  one  somewhat  akin  to  a 
forced  vortex;  what  is  equally  surprising  is  that  the 
total  temperature,  initially  uniform  at  the  inlet  and 
equal  to  97°F,  spontaneously  separated  into  hotter  stream 
of  about  118°F  near  the  outer  wall  and  colder  one  of  83°F 
near  the  inner  wall,  with  the  difference  as  distinct  as 
35°F!  This  latter  reminds  us  of  none  other  than  the  Ranque 
Hilsch  effect. 

Faced  with  the  severity  of  dynamic  flow  field, 
which  posed  a  serious  threat  to  the  subsequent  use  of  the 
rig  for  flutter  testing  of  airfoils,  the  annular  cascade 
was  modified  and  both  inner  and  outer  walls  were  provided 
with  tuned  acoustic  absorber  .  And  this  did  remove 
the  unacceptable  dynamic  flow  disturbances.  Ever  since, 
the  vehicle  has  successfully  been  in  use  for  aeroelastic 
purposes,  the  details  of  which  are,  however,  outside  of  our 
present  interest. 

Instead,  we  focus  our  capital  concern  to  what 
happens  to  the  profiles  of  steady  flow,  now  that  the  un¬ 
steady  fluctuation  has  been  eliminated.  The  answer:  the 
chance  in  the  velocity  and  temperature  distribution  has 


vanished.  The  free  vortex  remains  so  throughout  even  above 
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the  swirl,  where,  before  the  suppression  of  organized 
disturbance,  it  has  been  converted  into  a  forced  vortex 
type;  the  total  temperature  remains  uniform  in  the  radial 
direction  throughout  --  the  Ranque-Hilsch  effect  is  gone. 

Beyond  doubt  the  acoustic  streaming  did  somehow 
deform  the  steady  flow  field,  both  in  velocity  and  tem¬ 
perature,  and  this  affords  an  unmistakably  obvious  clue 
into  the  mechanism  of  little  understood  Ranque-  Hilsch 
effects,  which  we  shall  discuss  in  some  detail  below. 

We  recall  that  in  the  Ranque-Hilsch  tube  (Ranque 
1933;  Hilsch  1947),  the  compressed  air  enters  near  one 
end  of  a  single  straight  tube  through  one  or  several 
tangential  injection  nozzles.  Then, once  within  the  tube, 
the  swirling  air  segregates  by-  itself  into  two  streams  of 
different  total  temperature:  the  hotter  air  near  the 
periphery  of  the  tube  and  the  colder  one  at  the  centerline, 
a  separation  effect  already  mentioned  with  regard  to  the 
annular  cascade.  Between  the  Ranque-Hilsch  tube  and  the 
annular  cascade,  the  visible  difference  in  the  internal 
flow  passage  is  that  the  former  is  made  of  a  single  tube, 

*  We  eliminate  the  possibility  of  vortex  breakdown  (e.g. 

Hall  1972)  on  the  following  grounds.  First,  the  breakdown 
is  essentially  a  steady  phenomenon;  the  one  described  here 
was  unsteady.  Second,  the  measurement  in  the  annular  cascade 
did  not  exhibit  any  reversal  of  flow  in  the  axial  direction. 
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the  latter  of  an  annulus  (this  will  turn  out  to  be  a  not  so 
trifling  dissimilarity  as  it  might  seem  now).  As  a  matter 
of  further  detail,  in  the  conventional  Ranque-Hilsch  tube 
the  cold  air  is  immediately  extracted  from  an  orifice 
located  on  one  end, near  the  inlet  nozzle, and  the  hot  air 
spiralling  downstream  escapes  from  the  other  end  where  a 
throttling  exhaust  valve  is  located  \  this  is  the  so-called 
counter-flow  type.  Even  by  closing  the  cold  orifice,  the 
air  flowing  only  in  one  direction  toward  the  exhaust  valve 
can  still  produce  the  radial  separation;  this  is  called 
uni -flow  type. 

Detailed  measurements  of  the  internal  flow  distri¬ 
bution  in  the  Ranque-Hilsch  tube,  taken  at  the  condition  of 
optimam  cooling, show  that,  in  every  instance  a  forced  vor¬ 
tex  type  is  formed  immediately  near  the  entrance  to  the 
tube  ,  even  at  a  location  as  practically  close  as  possible 
to  the  inlet  nozzle  (for  uni -flow  type,  Eckert  and  Hartnett, 

1955,  Hartnett  and  Eckert  1957,  Lay  1959;  for  uni-flow  type 
with  vortex  chamber,  Savino  and  Ragsdale  1961;  for  counter¬ 
flow  type  Scheller  and  Brown  1957,  Sibulkin  1962,  Takahama 
1965,  Bruun  1969).  The  forced  vortex  occupies  the  almost 
entire  cross  section  (except  in  the  neighborhood  of  the 
boundary  layer  on  the  tube  periphery,  of  course)  and 
remains  so  from  the  entrance  to  the  exit.  Mark  with  atten¬ 
tion  that  at  this  condition  any  vestige  of  what  may  be  char¬ 
acterized  as  a  free  vortex  type  has  nowhere  been  detected.  Also  , 
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right  at  the  entry,  the  radial  separation  of  total 
temperature  occurs.  Contrary  to  some  earlier  belief,  the 
maximum  tangential  velocity  near  the  periphery  of  tubes 
needs  not  to  be  supersonic  to  create  the  effect,  even  the 
speed  of  500  ft/sec  or  so  suffices. 

Although  the  actual  total  temperature  separation 
in  Ranque-Hilsch  tubes  is  beyond  all  question,  none  of  the 
theoretical  explanations  devised  so  far  appear  to  have 
found  unreserved  acceptance.  Take,  for  example,  the  tur¬ 
bulent  migration  theory  (Van  Deemter  1952;  Deissler  and 
Perlmutter  1960;  Linders trdm-T  z  1971).  This  rests 
upon  the  assumption  that  when  a  lump  of  fluid  migrates 
radially  by  turbulent  motion,  it  tends  to  separate  the 
total  temperature  by  the  combination  of  the  following 
two  separate  mechanisms  of  stochastic  origin:  (1)  formation 
of  a  forced  vortex  and  (2)  creation  of  a  static  temperature  dis¬ 
tribution  approaching  an  adiabatic  one,  the  latter  through  the 
heat  transfer  process  in  a  centrifugal  field  originally 
postulated  by  Knoernschild  (1948).  However,  confrontation 
with  the  experimental  evidence  already  available  in  the 
literature  reveals,  that  the  contention  of  turbulence  as 
a  dominant  catalytic  agent  for  the  Ranque-Hilsch  effect  appears 
to  suffer  from  a  serious  flaw. 

Let  us  turn  our  attention  temporarily  away  from 
the  Ranque-Hilsch  tube  proper,  and  inspect  the  measurements 
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in  apparatus  with  the  tangential  injection  identical 
to  the  one  for  the  Ranque-Hilsch  tube,  but  constructed 
instead  to  create  a  vortex  for  purposes  other  than 
energy  separation  (Ter  Linden,  1949,  f«r  cyclone  separa¬ 
tor;  Keyes  1961,  for  containment  of  fission  material; 

Tsai  1964,  for  plasma  jet  generator;  Gyarmathy  1969, 
for  von  Ohain  swirl  chamber;  Batson  &Sforzini  1970,  for 
swirl  in  solid  propellant  rocket  motors.)  There,  unmis¬ 
takable  free  vortex  type  prevails  at  the  entrance  and 
elsewhere,  with  the  obvious  exception  of  the  innermost 
core  near  the  tube’s  centerline;  the  total  temperature 
at  the  entrance  remains  virtually  uniform  in  the  radial 
direction  and  equal  to  the  inlet  total  temperature  (Batson  & 

*To  be  more  precise ,  Batson&  Sforzini '  s  data  shows  that  the 
total  temperature  is  uniform  from  the  periphery  of  the 
tube  to  the  boundary  of  the  inner  core ,  which  occupies 
about  10%  of  tube  radius  from  the  centerline. (In  the  core, 
the  tangential  velocity  distribution  is  of  a  forced  vortex 
type  and  the  total  temperature  dips  slightly.)  Under 
conditions  fulfilled  in  their  experiments,  this  is  theoreti¬ 
cally  consistent  with  Mack's  results  (1960)  where  he  has 
shown  that  even  for  viscous,  heat-conducting  flow,  the  total 
temperature  of  a  free  vortex  remains  virtually  uniform  (and 
exactly  so  for  the  Prandtl  number  of  0.5)  provided  the 
swirl  Mach  number  is  less  than  one. 


Sforzini,  ibid.)  Even  in  these  test  rigs,  the  turbulence 
level  would  be  more  or  less  the  same  as  for  the  Ranque- 
Hilsch  tube.  Contrast  this  with  the  sudden  formation  of  a 
forced  vortex  and  the  separation  of  total  temperature  right 
at  the  entrance  of  the  Ranque-Hilsch  tube.  If  turbulence 
is  the  primary  agent  at  work,  then  under  circumstances  not 
unlike  each  other,  why,  in  the  particular  case  of  the 
Ranque-Hilsch  tube,  can  it  eradicate  any  traces  of  a  free 
vortex  and  suddenly  separate  total  temperature,  while  in 
the  others  it  can  'still  preserve  a  predominantly  free  vortex 
and  virtually  uniform  total  temperature? 

This  dichotonomous  branching  has  been  left  un¬ 
explained  by  the  turbulent  migration  theory.  Although  the 
space  does  not  permit  us  to  dwell  on  the  details  of  other 
theories  (Schepper  1951;  Sibulkin  ibid.),  they  do  not 
appear  to  be  free  of  similar  serious  objections. 

The  experimental  evidence  mentioned  in  the  open¬ 
ing  of  this  section  on  the  annular  cascade  compells  us  to 
turn  toward  acoustic  streaming  as  the  more  dominant  cause 
of  the  Ranque-Hilsch  effect  --  the  acoustic  streaming  in¬ 
duced  through  the  Reynolds  stresses  which  are  caused  by 
organized  periodic  disturbances  rather  than  by  stochastic 
motion . 

Close  scrutiny  of  the  available  past  literature  on 
the  Ranque-Hilsch  tube  reveals,  surely,  the  allusion  to  the 


presence  of  an  intense  periodic  disturbance  observed  by  many 
experimenters.  Hilsch(ibid. )  himself  mentions  that  a 
boiling  sound  was  audible  if  the  exhaust  valvewas  set  at 
optimum  position  for  cooling.  McGee  (1950),  Savino  and 
Ragsdale  (ibid.),  Ragsdale  (1961) .Kendall  (1962,  fora 
vortex  chamber),  and  Syred  and  Beer  (1972)  recount  in 
one  way  or  the  other,  the  disturbance  of  pure  tone  type, 
whistle  or  scream.  In  fact,  Savino  and  Ragsdale  record 
an  incident  where  a  loud  screaming  noise  was  accompanied  by 
10°  -  20°F  change  in  total  temperature,  a  phenomenon 
where  the  experience  of  the  annular  cascade  leaps  immediately 
to  mind.  None  of  them,  however,  proceeded  beyond  the  stage 
of  giving  passing  observations  to  it. 

To  a  certain  extent,  the  work  of  Sprenger  (1951) 
foreshadows  our  premises  in  its  spirit.  In  the  Ranque-Hilsch  tube 
with  its  hot  end  closed  and  only  its  cold  end  open,  he  measured 
periodic  disturbance  by  spreading  Lycopodium  to  form  a 
Kundt  pattern.'  However,  the  pattern  was  apparently  used  to 
measure  only  the  wave  length  of  discrete  disturbances,  since 
he  did  not  pin  the  Ranque-Hilsch  effect  down  to  the  acoustic 
streaming.  Rather,  by  appealing  to  the  analogy  of 

fC  ~ 

Reynolds  (1961) ,  while  advocating  the  turbulence  migration 
theory,  refers  Sprenger 's  idea  as  due  originally  to  Ackeret 
without  citing  the  reference:  to  date,  we  have  been  unable 
to  locate  the  original  source. 
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the  resonance  tube  (e.g.  Hartmann  1931),  he  later  simply 
suggested  (1954)  that  the  organized  unsteadiness  might 
produce  the  energy  separation. 

Highly  suggestive  also  were  the  circumstances 
which  led  to  the  discovery  of  vortex  whistle  by  Vonnegut 
(1954) .  While  engaged  in  experiments  exploring  the 
application  of  the  Ranque-Hilsch  cooling  effect  (as  a 
possible  means  of  measuring  the  true  static  temperature 
of  air  from  aircraft  in  flight) ,  Vonnegut  (1950)  observed 
the  presence  of  a  pure  tone  noise.  Although  he  did  not 
connect  it  with  a  mechanism  of  the  Ranque-Hilsch  effect, 
from  this  hint  he  constructed  a  musical  instrument,  the  so- 
called  vortex  whistle,  where  air,  injected  tangentially  into 
a  cylinder  of  larger  diameter,  swirls  into  a  smaller  tube; 
the  sound  thus  emitted  is  found  to  have  a  discrete  frequency, 
which  is  proportional  to  flow  rate.  Recall,  now,  that  the 
frequency  of  the  pure  tone  noise  in  the  similar,  swirling 
flow  within  the  annular  cascade  was  also  proportional  to  the 
Mach  number. 


Strickly  speaking,  finer  distinction  has  to  be  drawn  between 
the  two  frequency-swirl  relationships,  as  will  be  made  clear 
in  Section  7. 


1 . 2  Outline  of  Present  Investigations 

Against  the  precedent  setting,  we  shall,  in  the 
present  paper,  pose  the  following  model  problem:  periodic 
disturbances  in  swirling  flows  within  straight  co-annular 
cylinders  or  a  single  tube.  We  shall  solve  it  by  deriving 
an  explicit  expression  for  its  acoustic  streaming;  then, we 
shall  seek  to  display  such  key  features  as  the  transfigura¬ 
tion  of  steady  swirl  from  .one  type  to  another  at  certain 
threshold  steady  swirl;  and  we  shall  attempt  to  explain 
the  Ranque-Hilsch  effect  on  the  basis  of  streaming  caused 
by  Reynolds  stresses  due  to  organized  periodic  disturbances . 

Acoustic  streaming  is,  of  course,  an  induced 
steady  or  d.c.  component,  and  as  such  we  have  to  distinguish 
it  sharply  from  the  base,  steady  flow  initially  imposed  before 
the  disturbances  are  set  up.  For  brevity,  we  shall,  here 
and  henceforth,  refer  to  the  latter  simply  as  steady  flow 
and  its  sum  with  the  former  as  the  total  d.c.  component. 

Now,  without  a  single  exception,  the  only  known 
analytical  method  to  obtain  streaming,  the  present  one  not 
excepted,  is  to  resort  to  the  use  of  a  perturbation  scheme 
and  take  a  temporal  average  of  the  second-order  equation, 
which  contains  products  of  the  first  order  quantities.  Thus, 
if  we  started  from  the  conventional  boundary  layer  equations, 
which  corresponds  of  course  to  the  first  order  approximation 
to  the  full  Navier-Stokes  equations,  we  would  be  asked  as  to 
the  effects  of  what  are  collectively  called  the  higher-order 
approximation  to  boundary  layer  theory  (e.g.  Van  Dyke 
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1969)  on  the  streaming.  In  this  very  connection,  upon 
treating  the  problem  of  streaming  around  an  oscillating 
cylinder,  Stuart  (1966)  justly  voiced  a  note  of  caution 
on  the  possible  effect  of  curvature,  which  could  be  of 
the  same  second  order  as  the  streaming  itself.  (For 
this  particular  effect,  in  his  definitive  work  on  stream¬ 
ing  for  incompressible  flow  otherwise  in  a  state  of  rest, 
Riley  (1967)  shows  conclusively  from  a  matched  asymptotic 
expansion  that,  as  far  as  the  leading  term  of  the 
streaming  is  concerned,  the  curvature  has  no  influence 
within  the  unsteady  boundary  layer.)  In  the  present  case 
we  are  besieged  with  more  than  a  single  effect  of  possible 
second-order  correction.  For  example,  both  steady  and 
unsteady  boundary  layers  formed  over  the  cylindrical  sur¬ 
faces  present  the  problem  of  a  steady  as  well  as  an  unsteady 
displacement  thickness;  the  fluctuation  of  tenperature  gives 
rise  to  changes  in  the  viscosity,  which,  coupled  with  tem¬ 
poral  variation  in  strain,  might  beget  additional  Reynolds 
stresses,  as  will  be  found  to  be  indeed  the  case;  the  flow 
being  compressible,  even  the  effect  of  the  second  coefficient 
of  viscosity  must  be  assessed,  as  has,  in  fact,  been  done 
by  Van  Dyke  (1962-a)for  steady  compressible  boundary  layers  . 

To  face  these  problems  once  and  for  all,  we  shall 
abandon  the  standard  boundary  layer  equations  and  start 


afresh  with  the  full,  compressible  and  unsteady  Navier- 
Stokes  equations,  retaining  even  the  second  coefficient 
of  viscosity.  Under  the  conditions  of  several  parameters 
to  be  small,  we  shall  use  the  matched  asymptotic  expan¬ 
sions  to  ferret  out  the  leading  term  of  the  acoustic 
streaming  in  swirling  flow  within  co-annular  cylinders; 
by  following  this  avenue  of  plunging  into  the  equation 
in  its  full  generality,  we  can  not  avoid  somewhat  elabo¬ 
rate  alegbra,  which  constitutes  Section  3  through  6. 

One  of  our  centerpiece  results  is  equation  (42) , 
which  expresses  the  acoustic  streaming  in  the  circum¬ 
ferential  direction  near  the  cylindrical  surfaces .  This 
will  explicitly  show  the  following:  it  suffers  a  sudden 
reversal  of  its  direction  above  a  threshold  steady  swirl, 
the  physical  reason  being  due  to  the  Doppler  shift  caused 
by  swirl;  and  this  holds  regardless  of  the  values  of  the  stream¬ 
ing  Reynolds  number.  At  the  same  time,  the  absolute  mag¬ 
nitude  of  acoustic  streaming  itself  becomes  considerably 
increased. 

The  specific  value  of  swirl  at  the  threshold 
depends  on  the  explicit  relationship  between  the  frequency 
and  the  prescribed  radial  distribution  of  tangential  velocity. 


The  derivation  of  the  latter  relation  will  be  found  in 
its  entirety  in  the  Appendix.  With  the  aid  of  this, 
we  shall  be  in  a  position  to  discuss  the  behavior  around 
the  threshold  in  Section  7. 

For  the  steady  free  vortex  distribution  between  co- 
annular  cylinders,  which  corresponds  to  the  one  for  the 
annular  cascade,  the  threshold  steady  swirl  will  be 
shown  to  range  from  subsonic  to  supersonic  tangen¬ 
tial  Mach  number,  its  specific  values  being  strongly  dependent 
upon  the  ratio  of  outer  to  inner  radius  of  the  cylinder 
and  the  wave  modes.  Below  the  threshold,  the  tangential 
streaming  on  the  surface  of  the  inner  cylinder  is  in 
the  same  direction  as  the  steady  swirl;  however,  beyond 
this, the  streaming  abruptly  reverses  its  direction  and 
starts  to  retrogress  in  the  direction  opposite  to  steady 
swirl.  Surely,  then,  this  tends  to  decrease  the  total  d.c. 
component  of  circumferential  velocity;  its  reduction  being 
sizeable  near  the  threshold,  close  to  the  inner  cylinder 
the  radial  profile  is  converted  into  one  not  unlike  a 
forced  vortex.  This  behavior  appears  to  be  consistent 
with  the  observation  made  about  the  annular  cascade  in 
1.1. 

We  turn  now  to  the  steady  Rankine  vortex  dis¬ 
tribution  within  a  single  pipe,  which  corresponds  to  the 
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Ranque-Hilsch  tube.  There,  for  the  fundamental  mode 
of  disturbance,  any  amount  of  swirl  always  makes  the 
tangential  streaming  near  the  tube  periphery  rotate 
in  the  same  direction  as  the  steady  swirl  itself;  for 
such  a  wave,  no  threshold  swirl  exists.  Its  magnitude 
becoming  of  considerable  strength,  the  total  d.c.  component  of 
circumferential  velocity  in  a  free  vortex  region  is  increased  and  the 
entire  Rankine  vortex  is  now  converted  into  a  forced 
vortex; if  such  a  fundamental  mode  of  disturbance  is  not 
sufficiently  excited, the  Rankine  vortex  remains  vir¬ 
tually  unaffected.  (  In  all  these,  if  the  calculated 
streaming  becomes  of  considerable  magnitude,  obviously 
the  small  disturbance  approximations  break  down.  We 
assert,  however,  that  as  usual  this  does,  at  the  very 

least,  indicate  what  is  to  be  expected  in  real  situa- 

§ 

tions) .  In  the  light  of  our  discussion  in  1.1,  this  seems 
to  explain  the  appearance  of  swirl  either  in  the  form  of 
a  forced  vortex  preferred  for  the  Ranque-Hilsch  tube  or 
of  a  Rankine  vortex  for  others .  In  agreement  with  this 
also  is  the  following  observation  of  Takahama  and  Soga 
(1966)  made  in  a  Ranque-Hilsch  tube:  while,  at  the 
optimum  position  of  the  exhaust  valve,  the  forced  vortex 
is  found,  in  the  same  tube  with  the  exhaust  valve  now  at  off- 

optimum  position,  one  indeed  finds  a  well-defined  Rankine 


vortex. 


Whatever  the  initial  study  swirl  distribution 
may  be,  this  feature  they  both  have  in  common:  the  con¬ 
version  toward  forced  vortex  type  by  streaming  --  and 
this  tends  to  separate  the  flow,  with  the  initially  uniform 
total  temperature,  into  hotter  air  near  the  outer  radius  and 
colder  air  near  either  the  inner  radius  or  the  centerline; 
this  is  indeed  the  Ranque-Hilsch  effect. 

Details  of  the  physical  discussion  we  summarize 


in  Section  8. 


2.  Statement  of  Problem 


We  pose  the  problem  of  obtaining  the  acoustic 
streaming  in  an  annular  duct  between  two  circular  cylin¬ 
ders,  straight  and  stationery,  as  sketched  in  Figure  1:  vq 
denotes  the  radius  of  the  outer  cylinder,  r^  of  the  inner  cylinder. 
The  fluid  is  compressible  and  taken  as  a  perfect  gas  . 

We  assume  that,  throughout  the  entire  duct  length,  L, 
of  our  interest,  the  steady  boundary  layers  formed  over 
the  cylindrical  surfaces  are  thin;  in  the  inviscid  annular 
region  bounded  by  and  outside  of  them,  both  the  cir¬ 
cumferential  and  axial  velocity  are  the  predominant 
components  of  steady  flow,  as  shown  in  Figure  1.  Super¬ 
imposed  upon  this  steady  flow  are  the  fully  three-dimen¬ 
sional  unsteady  excitations,  whose  forms  are  to  be 
specified  in  the  inviscid  annular  region  and  whose 
streaming  effects  we  are  interested  in. 

Before  decomposing  into  steady  and  unsteady  parts, 
we  begin  by  writing  out  unsteady,  compressible  Navier- 
Stokes  equations  in  cylindrical  coordinates  (r,<j),z)  with 
corresponding  velocity  q  =  (u,  v,  w)  : 
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(2)  denotes  the  second  coefficient  of  viscosity,  the  one 
that  vanishes  for  a  monatomic  gas,  y  the  coefficient  of 
viscosity,  Pr  the  Prandtl  number,  0  temperature,  the  rest  of 
the  notation  being  standard.  The  inertia  terms  of  the  momentum 
equations  are  rendered  in  the  above  form  in  order  to  ease  the 
subsequent  streaming  calculations.  In  the  energy  equation, 
both  the  specific  heat,  c^,  and  Pr  are  taken,  as  usual,  to 
be  constant;  equation  (1-g)  simply  states  the  viscosity 
law. 

The  boundary  condition  on  the  cylindrical  surfaces  is 
q  =0.  The  thermal  condition  depends,  in  general,  on  the 
details  of  the  unsteady  heat  transfer  through  the  walls.  For 


21 


simplicity,  we  assume  that  either  the  wall  temperature 
is  maintained  to  be  constant,  its  unsteady  part  being 
kept  to  be  equal  to  zero,  or  the  walls  are  insulated. 
(Dunn  and  Lin  (1955)  have  shown  that  the  former  corres¬ 
ponds  to  the  situation  where  the  thermal  inertia  of  the 
wall  prevents  the  surface  temperature  from  responding 
to  high  frequency  fluctuation,  as  to  be  expected.) 


j 
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We  now  decompose  the  flow  field  into  the  steady  part  and 
the  one  related  to  unsteady  disturbance  such  as 


q  -  Q  +  q', 

where  Q  =  (U,  V,  W)  and  q'  ■ 
P  =  P  +  P'. 

..  o  =  R  +  p  "* , 
e  -  e  +  e 

y  =  M  +  \i' , 
y(2)=  M^2)  + 


(2-a) 

(u",  v',  w');.  likewise 

(2-b) 

(2-c) 

(2-d) 

(2-e) 

(2-f) 


Henceforth,  primes  denotes  unsteady  disturbances ,  vhich  contain  both 
a.c.  and  d.c.  components,  the  latter  being  streaming  in¬ 
duced  by  the  former. •  In  an  effort  to  seek  a  tractable 
ingress  to  the  problem,  we  asstime  disturbances  to  be 
organized,  regular  and  endowed  with  the  fundamental  frequency 
w,  the  stochastic  part  of  it  being  neglected  in  comparison, 
in  accordance  with  the  discussions  of  Section  1.1.  Of  our 
central  interest  is  the  effect  of  <§  upon  the  streaming. 

Henceforth,  it  is  convenient  to  render  r,  z,  and  t 
dimensionless,  by  referring  length  to  r^  and  time  to  «*>  ,  they 
being  denoted  as  starred  quantities;  i.e.  r*  =  r/r^,  z*  = 
z/r^,  and  t*  =  u>t.  We  also  define  the  ratio  of  the  outer 
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cylinder  radius  to  that  of  the  inner  cylinder 

^  «  X  .  (3) 

ri 

The  next  subsection  concerns  some  preliminary  con¬ 
sideration  on  the  steady  flow,  followed  by  the  one  on 
the  unsteady  disturbance. 
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3.1  Base,  Steady  Flow 

We  assume  that  the  steady  flow  is  axisymmetric 
and  apply  the  standard  outer  and  inner  expansion  (Van 
Dyke,  1962  a  and  b)  in  powers  of  R  where  Rg  is  the 

Reynolds  number,  whose  choice  of  characteristic  values 
is  left  unspecified  at  this  point.  We  first  examine 
the  leading  terms,  denoted  by  subscript  0, in  the  follow¬ 
ing  outer  expansions: 


■*  -1/2  ■* 

Q  -  Qo  +  Re  Ql  +— ' 

P  -  P0  +  Re"172?!  +---, 

R  -  R0  +  Re'1/2R1  +— • 

e  =  O0  +  Re_1/201+_““, 

M  =  H0  +  Re_1/2 

M<2>  =  +  Re-1/2M1(2^+— -, 


(4-a) 
(4-b) 
(4-c) 
(4-d) 
(4-e) 
(4-f ) 


where  Qq  =  (Uq,  Vq,  Wq)  ,  =  (U^,  Vj^,  W^) ,  etc.  As  usual, 

matching  yields  Uq  =  0  on  the  walls  and  henceforth  we  are 
interested  in  the  case  where 


u0  -  0,  (5) 

at  any  values  of  r*  standing  for  the  outer  variable.  All  the 
leading  terms  or  "inviscid"  flow  are  taken  to  be  independent 
of  z*(as  well  as  4>)  ,  in  view  of  the  assumption  of  thin 
boundary  layers  stated  in  Section  2;  they  are  all  functions 
of  r*only,  and  may  be  written  out  explicitly  as 
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-  Vr*>. 


W, 


=  W0(r*), 

=  Vr*>’ 


R, 


=  R0(r*), 


e. 


-  e0(r*), 


(5-b) 

(5-c) 

(5-d) 

(5-e) 

(5-f) 


etc.  Then,  the  only  nontrivial  equation  among  equation  (1) 
is  the  following, 
dP„ 


I  D  n  ^  _  y 

i*  0  0  3r*  , 


(6) 


which  expresses  the  radial  equilibrium  of  r-component  of 
momentum.  As  is  well  known  in  the  theory  of  turbomachines 
(e.g.  Marble  1964,  p.  144),  if  the  corresponding  stagnation 
enthalpy  be  uniform  in  the  radial  direction,  and  further¬ 
more  if  and  only  if  VgCr^Obe  of  free  vortex  type,  Wq  remains 
uniform  in  the  radial  direction;  for  other  steady  tangential 
velocity  distribution,  Wg  varies  in  the  radial  distribution 
even  for  uniform  stagnation  enthalpy. 

If  we  assume  that  the  corresponding  entropy  remains 
constant  everywhere , then 

Pn 


V 


=  constant, 


(7-a) 


or  equivalently 


86r 


>*0  ~TF*~  TF*  =  °. 


(7-b) 
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Finally  the  radial  dependence  of  the  inviscid  acoustic 

2 

speed,  defined  as  c  (y-  1)  0q  =  Aq,  is  given  by 

A 

[A0(r>v)]2  =  [Aq  (r*  =  X)]2-  (y-  l)J^[VQ)2  dr*,  (8) 

where  AQ(r*  =  X)  denotes  acoustic  speed  at  the  periphery 
of  the  outer  wall  in  the  co-annular  duct. 

The  leading  terms  of  the  inner  expansion,  con¬ 
stitutes,  of  course,  the  conventional,  compressible  steady 
boundary  layer  equations  (the  second  coefficient  of 
viscosity  does  not  show  itself  until  the  third  order 
(Van  Dyke  1962  a)).  On  both  cylindrical  surfaces,  the 
boundary  layers  are  rotationally  symmetric  and  develop 
in  the  z  direction.  Their  structures  being  complicated 
by  the  presence  of  two  components  of  inviscid  stream, 
swirl  and  axial  flow,  one  would  certainly  have  to  rest 
content  to  derive  them  by  such  methods  as  adopted,  for 
example,  by  Taylor  (1950)  for  a  swirl  atomizer  problem 
or  the  like.  However,  as  will  be  shown  later,  in  so  far 
as  we  limit  out  attention  to  the  leading  term  of  acoustic 
streaming  under  the  conditions  where  Rg  is  large 
and  the  steady  swirl  not  far  from  its  threshold  value, 
the  case  of  our  deepest  interest,  the  details  of  the 
steady  boundary  layer  need  not  be  worked  out. 
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3 . 2  Linearized,  Inviscid  Waves 

We  now  turn  to  q'  of  equation  (2-a)  and  the 
other  similar  unsteady  disturbances,  upon  which  our 
emphasis  lies . 

Before  ushering  in  the  more  precise  formalism 
of  the  asymptotic  expansion,  we  pave  the  way  for  it  by 
considering,  in  less  stringent  manner,  the  linearized 
inviscid  disturbances.  By  the  usual  small  perturbation 
of  equation  (1)  around  the  inviscid  steady  flows  of 
Section  3.1  and  tentatively  setting 

u'  -  u<r*)ei<"»'+  kriz*-  c*> 

v-  =  v(r*)ei<m*+  kriz*‘  c*> 

w'  =  kriz*‘  c*> 

and  the  like,  we  obtain  immediately 


(9-a) 

(9-b) 

(9-c) 


ifp  +  +  iktiRQW  -  0, 


(10-a) 


o  —  i  d-  —  dP  a 

ifu  -  V0v  =  -  ^  apr  +  ^2  dP^ 

ifv  +  ^  J^(r*v0)u  -  -  A  P. 


(10-b) 


(10-c) 
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0 

eo 


fr--  -1)  pE-  =  0 
Y  po 


(10-e) 


P 


(e0p 


+  Ro8), 


(10-f  ) 


with  the  boundary  condition 

u  =  0  at  r*=  1  and  r*=  X,  (10-g) 

where 

£ (r*)=  (-  »+  V0  ^  +  kW  \  r.  .  (10-h) 

In  the  above,  the  quantities  with  the  subscript  0  again 
stand  for  the  ones  of  steady  inviscid  flow,  all  being 
functions  of  r*only;  in  deriving  (10-b)  ,  the  radial 
equilibrium  relation  of  equation  (6)  has  been  taken  into 
account.  These  set  of  equations  are  idential  to  the  ones 
studies  by  Kerrebrock  (1977) . 

By  inspection  of  the  above,  we  realize  at  once 
that  the  real  parts  of  u'  and  of  all  the  rest  oscillate 
with  a  phase  difference  equal  to  --  a  fact  first  exploited 
by  Lord  Kelvin  (1880)  for  his  analysis  of  incompressible 
flow;  we  will  use  this  later  for  the  present  compressible 
flow  alike. 

Of  crucial  importance  among  the  foregoing  equa¬ 
tions  is  the  fact  that  frequency  w  does  not  appear  by 
itself,  but  arises  in  the  exclusive  form  of  (10-h).  That 
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is, 


f(r*)=  p  »  +  ^V0  (r*)+  k  WQ(r^  r. 

where  m  and  k  are  wave  numbers  in  the  circumferential 
and  axial  direction,  respectively;  obviously  this 
embodies  the  Doppler  shift  caused  by  the  steady  inviscid 
swirl  Vq,  and  the  axial  velocity  Wq  .  Notice  that  the 
shift  is  dependent  upon  the  radial  position  r*. 

In  order  to  obtain  the  acoustic  streaming, 
these  linear  disturbances  are,  of  course,  to  be  enforced, 
as  the  external  excitation  .upon  the  viscous  flow  near 
the  walls;  therefore,  in  forming  parameters  containing 
the  frequency,  the  relevant  frequency  to  be  used  there¬ 
in  is  not  bare  w,  but  the  above  f,  in  which  form  it  has 
recommended  itself. 

With  regard  to  the  above  set  of  equations 
themselves,  together  with  the  boundary  condition  equation 
(10-g) ,  they  determine  the  relationship  between  the 
frequency  parameter  f  and  the  steady  flow  field,  wave 
numbers  and  outer  and  inner  tube  radii;  upon  obtaining 
the  acoustic  streaming,  we  shall  later  need  such  explicit 
formulae.  For  now,  instead,  we  proceed  directly  to  obtain 
the  viscous  response  near  the  solid  walls,  relegating 
derivation  of  frequency  relationship  to  the  Appendix. 
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4.  Construction  of  Expansion  Series  and  Matching 

For  both  inner  and  outer  cylindrical  surfaces,  an 
analysis  can  be  carried  out  in  the  same  way.  To 
render  it  concrete,  we  choose  here,  as  henceforth, 
the  one  near  the  inner  cylinder  located  at  r*=  1, 
the  other  one  corresponding  to  the  outer  surface 
being  simply  obtained  by  the  obvious  replacement. 

Then, for  the  reason  just  stated  in  the  preceed- 
ing  section,  the  appropriate  frequency  parameter 
pertaining  to  such  a  study  is  (10-h)  evaluated  at 
r*=  1,  denoted  as  f^: 

fi  -  f  <r*-  »  *  (-"  +  Vex  +  k  Wex)  ri- 

where  subscript  "ex"  henceforth  stands  for  the  external 
flow  or  the  leading  term  of  the  outer  expansion  of  steady 
flow  evaluated  at  r*=  1:  namely, 

Vex  =  V**=  X>’  (12_a) 

Wex  =  ^(r*=  i),  (12-b) 

and  the  like;  we  note  that  they  are  all  independent  of 
z*as  well  as  <J>(and  of  r*,  of  course). 
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In  characterizing  the  unsteady  response  near  such 
a  solid  surface,  three .dimensionless  parameters  dis¬ 
tinguish  themselves.  The  first  is  defined  by 


(13) 


where  u  is  an  amplitude  parameter,  a  measure  of 
the  intensity  of  disturbances.  The  second,  6  ,  is  the 
ratio  of  the  unsteady  boundary  layer  thickness  6'  to 
r^,  where  6 '  is  defined  as 


S' 

Hence, 

6 


Mexri  1  ? 

Rex  |  f±  | 


(14) 

(15) 


The  absolute  form  of  f ^  is  needed  in  the  above ,  for  it  can  and 
does  switch  its  sign  depending  on  the  relative  magnitudes 
of  w  and  the  external  velocities .  If  V  =  W  =  0,  these 

CA  CX 

two  parameters  are  reduced  to  the  ones  employed  by  Riley 
(ibid.)  for  his  investigation  of  streaming  in  the  absence 
of  base,  steady  flow;  there,  as  long  as  the  frequency  is 
prescribed  and  held  fixed,  both  remain  obviously  unchanged. 
This  is  to  be  contrasted  with  the  following  present  situation 
as  steady  flow  varies,  the  unsteady  boundary  thickness  6' 
changes;  consequently, 3  does  not  remain  the  same  even  for 


a  given  frequency.  For  example, 6'  or  3  increases  when 
is  decreased  by  the  change  in  Vex  only.  Likewise  a 
does  not  remain  the  same . 


The  third  and  last  parameter  e  is  a  measure  of 
relative  magnitude  between  steady  boundary  layer  thick¬ 
ness,  6  to  the  unsteady  one,  6';  the  former  is  defined 


(16-a) 


in  terms  of  the  duct  length,  L,  and  where  Re>  the  Reynolds 
number,  so  far  purposely  left  undefined,  now  denotes 
precisely 


Wexpex 

Re  = 

ex 


Hence , 


(16-b) 


(16-c) 


on  account  of  equations  (16-a)  and  (14) .  Recognize  here 
that  if  the  Reynolds  number  is  sufficiently  high,  e  can 
indeed  by  a  small  number,  in  so  far  as  we  are  interested  • 
in  the  region  not  too  remote  from  the  turning  point  to  be 
defined  later  as  f.  =0  . 


33 


Our  interest  is  directed  to  the  situation  where  a  , 
g  and  t  are  all  small  in  their  magnitudes.  Physically 
this  may  be  considered,  as  follows:  the  amplitude  of 
unsteady  excitation  is  small;  and  compared  to  the 
cylinder  radius,  the  unsteady  boundary  layer  thickness 
is  thin,  within  which  thinner  steady  boundary  layers  at 
the  high  Reynolds  number  is  embedded. 

We  name  the  annular  region  confined  between  the 
outer  edge  of  the  unsteady  boundary  layer  and  that 

'fc 

of  the  steady  boundary  layer  as  the  "middle  deck"; 
then,  the  "lower  deck"  naturally  suggests  itself  as  the 
annular  cross  section  between  the  outer  edge  of  the  steady 
boundary  layer  and  the  cylindrical  solid  surface.  The 


The  main  stream  consisting  of  both  swirl  and  axial 
velocity,  to  be  precise  we  have  to  take  the  possible  pre¬ 
sence  of  more  than  a  single  thickness  of  the  steady  boundary 
layer  into  consideration  (Cooke  1952) .  If  definiteness 
requires  it,  we  may  choose  the  thickest  among  them  as 
corresponding  to  the  outer  edge  of  the  steady  boundary 
layer. 
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double  deck  structures  are  present  near  both  the  outer 
and  inner  cylindrical  surfaces;  Figure  2  shows  this 
schematically  and  in  disproportionately  magnified  manner 
The  annulus  bounded  by  two  middle  decks  we  call  "core" 
We  turn  now  to  the  construction  of  series  expansions 
appropriate  for  each  region,  starting  from  the  middle 
deck. 
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4.1  Middle-Deck  Series  Expansion 

Let  n  be  a  middle  deck  variable  or  inner 
variable,  scaled  to  the  unsteady  boundary  layer  thick¬ 
ness  and  defined  by 


6 

We  first  turn  to  the  steady  flow  field  within 
the  middle  deck,  corresponding  to  the  first  terms  on  the 
right  hand  side  of  equation  (2) ,  and  express  them  in  the 


following  form: 

U  =  Wexn(r*,z*;  Re)  ,  (18-a) 

V  -  Vex  [1  +  e(r*.  z*;  Re)]  (18-b) 

W  =  Wex  [1  +  c  (r*,z*;  Re],  (18-c) 

R  =  Rex  X(t*,z*;Re)  .  (18-d) 

0  "  0ex  x(r*>  z*:  Re*  (18-e) 

P  =  RexVex2  ’H1**  z*JRe>  •  (18"f) 

M  =  Mex  a(r*,z*;  Re),  (18-g) 

and  M(2)  =  Mexa  z*;Re).  (18-h) 


■  -  -  -  ■_ 


By  definition,  the  middle  deck  lies  exterior  to  the 

steady  boundary  layer;  hence,  such  outer  expansions 

of  steady  flow  as  equation  (4)  and  (5)  are  applicable 

-1/2 

here.  We  expand  £,n  .  etc.,  in  powers  of  Rg  , 
which  is  proportional  to  c  $  from  equation  (16-c)  and 
obtain 


u  =  Wexe6n1(^v,z^!)+ -  * 

(19-a) 

V  =  w 

ex 

1  +  £Q(r*)+  e$51(i*,z*)+--- 1  > 

(19-b) 

W  =  W 

ex 

1  +  £0(r*)+  eB^O^.z*)*---] 

(19-c) 

R  -  R 

ex 

[  X0(^)  +  eBx1(^Z*>+--“] 

(19-d) 

e  -  e 

ex  | 

[t0(r*)  +  £BT1(r*fz*>4— -]  , 

(19-e) 

etc . .where 

£0  (r*  *  1)  =  CQ(r*=  1)  =0  ;  XQ<r*=  1)  =  t0(r*«  1)  --  *1. 

(19-f) 

The  last  one,  (19-f ),  follows  owing  to  the  very  definition 
of  quantities  with  subscript  "ex".  We  .then .expand  the  lead¬ 
ing  term  with  subscript  "0"  around  r*=  1,  make  use  of 
equation  (19-f)  and  finally  express  these  in  terms  of  the 
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middle-deck  variable  n.  Thus  we  have 


U  =  Wexe6  n1<n.z*>+---  . 

V  «  Vex  £  1  +  ajng  4-  a2(ng)2  + - 

W  =  Wex  [l  +  b^B  +  b2(ng)2  +  — 

R  =  Rex  [l  +  c^ng  +  c2(ng)2  + - 

0  =  0ex  [l  +  d-^ng  +  d2(ng)2  + - 

and  the  like,  where  a-^,  a2,  b^, 

constants . 


(20-a) 

+  £g?1(n,#)+ - J,  (20-  b) 

+  eg^Cn.z*)  + - J  ,  (20-c) 

+  egx1(n,z^)+ - ],  (20-d) 

+  ega^ (n, z*)  4 - J,  (20-e) 

b2,  ...,etc.  are  all 


For  the  unsteady  parts,  denoted  with  primes  in 
equation  (2) ,  we  expand  them  in  terms  of  a, g  and  e  with  the 
triple  indices  affixed  to  each  term:  for  instance, 


where  i  stands  for  the  order  of  a,  j  for  g,  and  e  for  k. 
Accordingly  we  write 


U'  -  ug  <U0'00  +  <XU'00  +  gu510  +  eUQ01  +—  )  , 
v"  =  u  ^v000  +  aV100  +  &vnm  +  evnm  4  )  * 


010 


001 


(21-a) 


(21-b) 


=  u  <w000  +  aw100  +  ewnm  +  EWnm  +  . 


010 


'001 


(21-c) 


P 

0  •* 


U  Rex  f i  (p'000+ap100  +  pp010  +  ep001  +"""  > 
-  Rex 

u  A  (p000+  ap100  +  Bp010  +  ep001  +""_  ^  • 


U  (0OOO  +  aGioo  +  B0oio  +  e0ooi  +—  >* 

CA 


(21-d) 

(21-e) 

(21-f) 
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*  u  A ~  000  +  ay  100  +  g  M  010  +  e  y  001  +  ^ 

ex 

u(2).  _  :  jex  ,  ,  (2)  >  (2)  .  (2)  .  _  u ^  (2). 

M  “  u  A  ^  000  +  ay  100  +  p  y  010  +  e  y  001  + 

ex 


w  +-), 
(21-h) 


where  in  the  last  two  equations,  the  leading  terms  are 
related  to  the  one  of  6'  on  account  of  (1-g)  as 


* 

do 

000 

dT 

W2* 

000 

dx 

(22-a) 

(22-b) 


The  derivatives  of  a  and  o^  '  defined  in  equations  (18-g) 
and  (18-h) ,  respectively,  are  to  be  evaluated  at  the  steady 
temperature  corresponding  to  the  inviscid  value  on  the  wall. 
Note  that  the  leading  term  of  u'  is  0(8) ;  the  others  are  taken 
to  be  0(1).  All  the  first  terras  in  the  brackets  with  indices 
(0,0,0)  representing  linearized  viscous  disturbances,  we  ex¬ 
press  them  explicitly  as 


u'  =  77  M  +  kr..z*-  t*) 

u000  u000  '•n;  e  1  * 

v  *  =  v  fnl  +  kr  .  z*-  t*) 

000  000  e  1  . 

w  •*  =  w  +  kr. z*-  t*) 

000  w000  '■n;  e  1 

and  the  like. 


(23-a) 

(23-b) 

(23-c) 


39 


4 . 2  Core  Series  Expansion 

In  the  core,  i*  ,  is  the  appropriate  radial 
coordinate  or  the  outer  variable.  The  expression  for 
the  steady  flow,  equation  (20) ,  is  immediately  applica¬ 
ble  to  the  core  as  well. 

For  the  unsteady  part,  we  expand  them  in 
series  structually  similar  to  (21)  except  for  u'  where 
the  leading  term  is  now  elevated  from  0(3)  to  0(1),  and 
equal  to  the  others .  Thus  we  have 


'000 


-  R 

I 

U  A 

~  e 

u  s 

~  M 


ex 

ex 

ex 

ex 

ex 

L 

ex 

M 


+  aU{0Q 

+  BUoio 

+  eUooi 

+--)  , 

(24-a) 

+  aV100 

+  evoio 

+  eVooi 

+ — )  , 

(24-b) 

+  aW£00 

+  eWoio 

+  eWooi 

+ — )  , 

(24-c) 

(Fooo  + 

aP100  + 

ePoio  + 

ep6oi+-->. 

(24-d) 

(Rooo  + 

aR100  + 

BR010  + 

eR001+-""^  * 

(24-e) 

(0ooo  + 

a0ioo  + 

B0oio  + 

c9001+~")  * 

(24-f ) 

<M6oO  + 

aM100  + 

BM010  + 

eM6oi+--> * 

(24-g) 

=  u  ^  <m6oo(2>  +  aM^™(2)  +  ^-(2)  +  cm;  (2) 


ex 


100 


mo 


‘001 
(24-h) 


where  again  the  leading  terms  of  the  last  two  equations  are 
related  to  ©000’  t'*ie  ^^rst  term  of  O'",  by 
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M000  "  GOOo[di  ]T= 


M 


000 


(2)  , 


0(r*>> 
(2)' 


=  e  -  fda^n 

y000  dx 

L  -I  T=  T 


,(r*)  . 


(25-a) 

(25-b) 


The  derivatives  are  to  be  evaluated  at  the  steady  temperature 
corresponding  to  the  inviscid  value. 

The  leading  terms  represent  the  linearized  inviscid 
disturbances  and  take  the  following  forms : 


000 

=  U  (**)  e"^ 

u000^  ; 

+  kr^z* 

-  t*), 

(26-a) 

000 

-  ''000<r*)e1<I"f 

+  kr . z* 

l 

-  t*), 

(26-b) 

000 

■  Wooo(r*>ei(m* 

+  kr±z* 

-  t*), 

(26-c) 

etc;  these  correspond  to  equation  (9)  but,  here  and  hence¬ 
forth,  we  are  denoting  them  by  the  formalized  triple  indices 
instead. 
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4 . 3  Matching  between  and  Middle  Deck  and  the  Core 

Standard  matching  procedures  (Van  Dyke,  1962 
a  and  b)  applied  between  the  middle  deck  and  the  core 
expansion  series  yields  at  once  the  required  matching 


condition.  For  u',  we  have 

Up00(l,<|>,z*,t*)=  0  ,  (27 -a) 

U001(1,l|),Z*,t:*)=  0  (27 -b) 

U'00(l,<j,,z*,t*)=  0,  (27 -c) 

UQlo(l,<0,z*,t*)=  lim  [uQOO(n,0>,z*,t*)-^^  n]  .  (27 -d) 

n  ->  oo 


The  last  equation  represents  the  influence  of  an. unsteady 
displacement  thickness,  exactly  analogous  to  its  steady 
counterpart  (Van  Dyke,  ibid.).  For  v",  we  have 


vooo 

(1 ,4>  , Z* , t*)  =  lim  VQ00(n,<})  ,z*,t*) , 

(28 -a) 

n  “ 

V001 

(1,4>  .  z*,t*)=  lim  VQ01(n,  4>,z*,t*)  , 

(28 -b) 

n  *■  oo 

V100 

(1 ,4>  ,  z*,t*)  =  lim  vC  00  (n ,  4>,  z* ,  t*) , 

(28 -c) 

n  -*■  oo 

V010 

(1,4>  ,z*,t*)=  lim  vQio(n,  - 

n  ( 

Koo 

H 
C  II 

t 

(28  -d) 

3r^ 

The  identically  same  relationships  hold  for  w',  p'.p'',  and 
0  "* .  The  matching  for  \i'  yields 
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M600(l.*»z*'t*)  =lim  . 

(29-a)  J: 

n-+°° 

i 

M6oi*1,*,z*,t1^=lim  ^mCn^.^.t*)  • 

n->co 

/ 

(29-b) 

M510(l,<|),z*>t*)*lim  |  yQ10(n,<|),z*,t*) 

-(Sj,  p] 

M-  I  ■ 

(29-c)  ! 

(2) 

TVia  Rama  mat'r'hino  rnnHiM'rmc  hnlrl  frrr  n  '  '  ' 

• 

The  same  matching  conditions  hold  for  y 
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4 . 4  Lower  Deck  and  the  Transfer  of  Boundary 

Conditions 

Even  in  the  lower  deck,  the  appropriate  radial 
variable  is  "n"  introduced  in  equation  (17)  .  Compli¬ 
cations  arise,  however,  in  the  lower  deck  due  to  the 
steady  flow  field  immersed  in  its  entirety  within  the 
steady  boundary  layer,  whose  structural  complexity  has 
already  been  mentioned  in  Section  3.1.  Though  rotational- 
ly  symmetric,  not  only  does  the  steady  profile  change 
radially  from  the  solid  wall  to  the  outer  edge  of  the 
steady  boundary  layer  but  also,  owing  to  its  axial 
development,  it  is  dependent  on  the  z  coordinate  as 
well;  due  to  the  presence  of  both  swirl  and  axial  flow 
in  the  main  stream,  the  steady  flow  field  is  fully  three- 
dimensional;  it  is  through  this  that  the  unsteady  dis¬ 
turbances  propagate. 

One  can  circumvent  this  difficulty  associated 
with  the  lower  deck  by  transferring  the  boundary  con¬ 
dition  on  the  solid  wall  to  the  inner  edge  of  the 
middle  deck;  such  bypassing  is  valid  for  terms  like 
(0,0,0),  (0,1,0)  and  (1,0,0)  but  not  for  (0,0,1).  Take, 

^(e) 

for  example,  v'  .  We  compare  its  exact  solution  denoted  as  v  , 
assumed  to  be  somehow  kn’own  in  both  middle  and  lower 
decks,  with  the  present  one  for  the  middle  deck,  denoted  here 
temporarily  as  Then  at  the  inner  edge  of  the 

middle  deck  located  at  n  =  e(z*) ,we  have 
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v'(e)(n  «  e  (z*),  z%,t*)«  v'(md)  (n  =  e(z*)  ,  z*,  <j>,t*) 


Substitute  equation  (21-b)  into  the  right  hand  side 
and  expand  both  sides  around  n  =  0.  Then  from  the 
no  slip  condition  on  the  wall,  the  left  hand  side 
becomes  0(e)  and  the  right  hand  side  yields 

v000  =  v010  =  v100  0  at  n  =  0. 


(30) 


Although  not  so  with  VqQ1  ,  this  will  not  present  diffi¬ 
culty  as  far  as  the  leading  term  of  the  streaming  is  con¬ 
cerned.  Exactly  identical  relationships  can  be  derived 

/p\ 

for  u^  and  w^(For  the  former,  match  3u'v  '  with  its  middle 
deck  counterpart,  equation  (21-a)  .)  Likewise ,  for  the  thermal 
condition  where  the  unsteady  temperature  on  the  wall  is 
maintained  to  be  zero,  wc  have 

6000  =  e010  =  9100  =  0  at  n  =  0,  (31) 


but 


Substitution  of  the  expansion  series  in  the  preceeding 
Section  4  into  equation  (1)  yields  sets  of  equations 
appropriate  for  each  order.  In  this  section,  we  assemble 
them  in  hierarchic  order,  and  describe  their  salient  features 


before  presenting  the  details  of  the  solution. 
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5.1  (0,0,0)  Core 

The  set  of  equations  is  given  by: 


!V  T  *  ^  1  kr.W  1 

1 +m-fr  L5°  '  ^  (1  +  co)J  +"TTL  ^ojiRooo 

+  ~^i  |  **  CrX°  ^000]  +  ir*x0V000  +  lkrixOWOOoj  =0, 


■-  (i  +  c0)  uooo  +  TT^  couooo 


0  |iUooo  +  im  if  [co  -  U  +  *o>]  u 

)  V  2  1 

f  ex _  1 

)"  fiAex  ? 


ikr.W 


(33-a) 


ex  1.  /I  J.  £  \ 

T~  r*  u  +  0'  v000( 


r*  R000 


dP 


000 

dr* 


(1  +  Cn>2 


(33-b) 


I* 


xo  f  ivooo  +  im  ~T7 

v 


“It  ‘  if*  (1  +  V  vooo 


+  —4W-e-  Cnv 


77 


0V000 


x  «  ^0+Vex  1+50  -  1_  im 

“  u000  3r*+  f±  ~ r*~  u000  (  r 

{V  f  *,  1  -  ikr.W 

iW000  +  im  ~T7  [^0  "  ^  (1  +  VJ  W000  +  -f-  ^0W 


r*  *000, 
(33-c) 
000 


W  3  Cn  -  ) 

ex  IT  l  =  _ikr 

T7  ^7*  u000  (  lkri 


000 


f . 

X0  e000  "  A”  P000  =  0  » 


(33-d) 

(33-e) 


“  Y  p000  =  t0R000  +  fy9  000’ 


ex 


(33-f) 
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with  the  boundary  condition 

U000(r=l)  =  0,  (33-g) 

and  the  similar  condition  to  be  satisfied  on  the  outer 
cylinder . 

This  set  of  equations,  once  restored  in  the 
original  variables,  is  identical  to  equation  (10) 
corresponding  to  the  linearized,  inviscid  disturbance. 

They  are  complete,  independent  of  (0,0,0)  of  the  middle 
deck  or  any  other  higher  order  equations.  They  being 
treated  in  detail  in  the  Appendix,  we  do  not  deal  with 
them  any  more  here. 
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5.2 

(0,0,0)  Middle  Deck 

The  set  of  equations  is  given  by: 

ip  + 

000 

Aex  |~du000  ,  .  -  ,  -  1 

ft  [  dn  +  lmv000  +  lkriw000j  '  °- 

(34-a) 

iv000  = 

4v  .  d2v000 

lV000  (r  1)1  ~^2 

(34-b) 

iwooo  = 

ikr.  _  *  d^w  . 

.  vooo<r=1>  1  -Jr  ■ 

(34-c) 

10OOO  + 

1  <T-l)  v“ooo  -  1  p1  — T  • 

(34-d) 

fi 

Y  V000(r=1)  =  p000  +  ®000  ’ 

(34-e) 

mAex 

where  +  corresponds  to  L  >  0  or  <  0,  respectively. 
The  boundary  conditions  are: 


=  v  =  w 

000  000 

-  0  , 

at  n  =  0 

(34-f) 

-  0  or  d^000 

°-  or  dn 

■  -  0  , 

at  n  =  0 

(34-g) 

o 
o 
o 
1  > 

i: 

8 

II 

C 

s/ 

(r*l)  , 

(34-h) 

kr-  -  * 

*>  =  -E  vooo(r=1)  • 

inA^  (Y'1)  Wr=1>’ 


m^xV000(r_1) 

m  V000(r=1)  ’ 


(34-i) 
(34- j ) 
(34-k) 
(34-1) 


-iff, 
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In  the  above,  the  radial  component  of  the  equation  of 
motion 

dPooo  _  n 

“  dn  u’ 

has  been  spent,  upon  combining  with  (0,0,0) order  of  core 
equations,  to  replace  the  pressure  terms  on  the  right 
hand  with  the  inviscid  tangential  velocity  on  the  wall, 
V000Cr=l).  Once  VqqqC^I)  is  specified,  the  set  is 
complete  and  can  be  determined  independently  of  any  other 
orders .  Absent  from  them  are  the  influences  of  varying 
steady  flow  field,  curvature,  the  second  coefficient  of 
viscosity  or  any  other  higher-order  effect  to  boundary 
layer  theory.  In  fact,  they  are  nothing  but  those  of  a 
conventional  compressible,  unsteady  boundary  layer, 
rendered  into  linearized  form.  Convective  terms  being 
eliminated,  the  equations  are  described  relative  to  a 
helically  moving  frame  of  reference:  a  coordinate  in 
screw  motion,  revolving  with  the  inviscid  swirl  Vgx  and, 
at  the  same  time,  advancing  with  the  inviscid  axial  velocity 
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5 . 3  (1,0,0)  Middle  Deck  --  Streaming  Part 

By  taking  the  temporal  average  of  (1,0,0) 
order  of  the  middle  deck, we  obtain  the  equations  which 
will  turn  out  to  be  the  leading  terms  of  streaming. 

We  first  list  the  ones  corresponding  to  continuity, 
tangential  and  axial  momentum  and  energy,  setting  aside, 
for  the  time  being,  radial  momentum  and  equations  of 
state: 


ex 


a  <u£oo^ 


3n 


+  3n  <p000  *  “OOO^  °’ 


(35-a) 


3n  <u000 


V6o°  <Vl00>+'^  3k(^600-3^>)]  , 


(35-b) 


&  <u6oo  •  w6oo>=±  (^r  <wi00>+ A  $v'J%SL>) 

BX 


(35-c) 


9  6  000 


<u000  3n  >  +  <v000  3^  (9000  "  (y"1)  Ae^  POOo)> 


+  ^000  3z*  (0ooo 

f . 


(y-!)  X"~  Pooo^ 


ex 


i  <P6oo  •  Real  [iT000.l«"*  + 

1  r  82  -  v  4.  fi  3  ( /  u  '  36  000  N 

pI  rr  <Gioo>  +  r~  3?r  v*  ooo  ~w~> 

r  L  3n  ex 


i 


1  (Y-l>  X — 

ex 


i  [<  (-^°)2>  *  . 


(35 -d) 


where  <  >  denotes  the  temporal  average  and  +  corresponds 
to  either  >  0  or  <  0.  The  boundary  .conditions  are 
as  follows 


<u100'>  “  ^v100^  ^w100^  “  0  at  n  -  0, 


(35-e) 


8  0”" 

<G1'oo>  -  0  or  -  0  at  n  -  o, 


(35-f ) 


and  for  n-+» 


Um  9<“l'0°> 

n-><x»  8n  n-voo  an 


8^e100^ 

an 


(35-g) 


As  for  the  boundary  condition  for  n^»  ,  the 
original  matching  condition  for  vfoo  rea<^s 

<v100>  =  <V100>r*l 

n-x» 

and  streaming  in  the  core  will  be  in  need  of  this;  for  now, 

as  usual  (Riley,  ibid.;  Batchelor  1967,  p.360.)  the  weaker 

3<v'l00> 

condition  <v{oo^n-*<>°  =  or  equivalently  — ^ —  "*■  0 

as  n->oo  suffices;  the  similar  conditions  apply  to  w£qq  and  6{qq; 
as  for  Ujqq,  one  does  not  need  any  condition  at  n  +  », 


<u£oo> ’ <v{oo>  ’  <w100>  and  <9100>are  comPletely 

and  uniquely  determined  in  the  middle  deck  from  the  above. 

These  d.c.  components  of  velocity  and  temperature  field  are 

dependent  upon  only  (0,0,0)  order  in  the  middle  deck  and 

they  are  independent,  among  others,  of  the  streaming  part 

for  (1,0,0)  in  the  core;  on  the  other  hand,  the  latter  will 

be  dependent  upon  the  former.  This  holds  regardless  of  the 

magnitude  of  the  streaming  Reynolds  number,  R  . 
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The  sequential  dependence  on  these  is  depicted 
schematically  in  Figure  3  as  a  "family  tree"  of  stream¬ 
ing  (applicable  to  d.c.  component  of  temperature  as  well). 

The  set  of  the  foregoing  equations,  (35),  is  independent  of 
variations  in  the  steady  state  and  of  both  steady  and  unsteady  dis¬ 
placement  effects.  They  are  independent  of  curvature 
effect,  as  can  be  readily  recognized  by  re-expressing 
them  first  in  terms  of  the  conventional  boundary  layer 
coordinates  and  identifying  them  as  time-averaged  forms 
of  the  boundary  layer  equations ,  where  the  second-order 
terms  of  the  unsteady  disturbances  are  retained:  The 
issue  of  curvature  effect  was  first  raised  by  Stuart  (ibid.) 
who  also  justified  its  neglect,  irrespective  of  the  value  of 
R  ,  on  physical  grounds  and  Riley  (ibid.)  confirmed  this 
afresh  by  a  matched  asymptotic  expansion. 

Observe  in  (35-b)  and  (35-c)  the  presence  of 

^v000 

terms  in  the  form  of  I-^qq  — ,  which  correspond  to 
the  Reynolds  stresses,  generated  by  the  fluctuation  of 
viscosity  responding  to  the  unsteadiness  in  the  tempera¬ 


ture  p 


and  coupled  with  the  temporal  oscillation  of  the 


strain  v000  ;  likewise  in  (35-d)  the  term  p' 


represents 


the  heat  flux  caused  by  a  similar  fluctuation  of  thermal 
conductivity  and  temperature  gradient. 

In  contrast  with  the  d.c.  component  of  velocity 
and  temperature,  the  temporal  averages  of  pressure  and 


density,  <p^QQ>and  <p{qo>>  respectively,  can  not  be 
determined  from  the  following  equations  and  boundary 
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conditions  alone,  which  have  been  intentionally  held  over 
so  far: 


9n  <Pioo>  0  * 


t . 

7  7T  <pioo>  =  <pfoo>  +  <eioo>  +  A~  <p6oo6ooo> 

ex  ex 

with 


(3  5-h) 

(35-i) 


lim  <P100>  =  <R100>rll  and  lim<p100>  =  <p{00>r*l. 
n->oo  n-*-°° 

They  need  to  be  matched  with  the  corresponding  ones  in  the 


core . 
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5 . 4  (1,0,0)  Core  -  Streaming  Part 

In  the  core,  (1,0,0)  streaming  is  driven  by 
the  streaming  of  the  same  order  in  the  middle  deck-- 
regardless  of  the  value  of  Rg ,  no  matter  how  large  or  small.  As 
is  now  well  know,  the  structure  of  streaming  outside  of  the  un¬ 
steady  boundary  layer  is,  however,  strongly  dependent 
upon  R  ,  as  pointed  out  first  by  Stuart  (1963,  p.  384) 

o 

and  subsequently  confirmed  explicitly  by  such  experiments 

as  for  a  cylinder  oscillating  in  otherwise  stationary 

fluid  (Bertelsen,  1974) . 

If  we  take  simply  the  temporal  average  of 

equations  of  the  order  of  (1,0,0)  in  the  core,  they 

correspond  to  the  case  of  R  <<  1  and  are  given  by: 

s 


<u100>  ■  °' 


*0  T7  I*  <!  +  «0>  <vi00> 


(36-a) 


vex2  1  2 

La  +  ?n><Ri"nn> 


f  .A  r 
r  ex 


+  r*  8r*[r*x0  <U0002>] 

r*x0<'V000  •>"  2  A~^  ^  +  V  ^000  'V600> 

ex 

=  _  _JL  Sp'  \ 

3r*  '*100', 


(36-b) 
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f . 


f . 


T“  Y  <P100>  T0  <R100>  +  X0  <e100>  +  A  1  <R000  *  6000  ^  ' 

ex  ex 


with  boundary  conditions 


(36-c) 


<V100>  r*=l  llm  <v100'> 

n->-co 


(36-d) 


^100^  r*=l~  Pim  > 


n-*-°° 


100' 


<01 00>  r*=l  llra  <5100> 

n->« 


(36-e) 

(36-f ) 


together  with  the  similar  boundary  conditions  to  be 
satisfied  on  the  outer  cylinder  surface.  In  deriving 
these,  we  have  made  use  of  the  phase  properties  of  (0,0,0) 
order  of  the  core  flow,  already  referred  to  in  Section  3.Z. 

Missing  from  the  above  set  is  the  energy  equation, 
which  vanishes  identically  by  virtue  of  equations  (7-b) , 
(33.-e)  and  (33-f)  :  the  equation  of  continuity  and  the 
momentum  equations  in  <f>  and  z  directions,  they  all  collapse 
to  yield  (36-a) . 

Viscous  terms  being  also  absent  in  the  above  momen¬ 
tum  equations,  they  are  of  inviscid  type.  The  boundary 
conditions  (36-d)  and  (36-e)  correspond,  however,  to  the 
non-slip  ones  on  the  "walls:  they  are  set  in  spiralling 
motion  by  streaming  in  the  middle  deck  and  to  the 
prescribed  surface  temperature  of  (36  -f)  •  Then,  not  only 
some  of  the  equations  are  lacking,  but  due  to  their 
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depressed  order,  ostensibly  they  do  uot  satisfy  the 
boundary  conditions.  The  latter  is,  however,  not  an 
indication  of  some  kind  of  breakdown,  but  the  missing 
ones  are  to  be  supplied  from  the  terms  of  higher  order 
than  the  present  (1,0,0).  That  this  is  so  can  be  illu¬ 
strated  by  referring  to  the  streaming  induced  by  fluc¬ 
tuation  in  incompressible  flow,  which  is  otherwise  at 

rest.  For  V  =0,  (36-b)  is  reduced  to 
ex 


1 

r* 


a 

ar* 


[r*x0<u6oo  >] 


'  i-«o<vooo2>‘- 


a 

3r* 


<  pioo>  ,  (37 ) 


and  we  note  that  <V£QQ>or  the  tangential  velocity  is 
now  completely  lost  from  this  or,  for  that  matter,  from 
everywhere  else.  It  can  be  recovered  only  from  the  con¬ 
sideration  of  viscous  terms  of  (1,2,0)  order,  a  higher  order 
term  in  the  expansion  series:  after  some  algebra,  we 
obtain 


a 

a?* 


a 

ar 


R<W1  <V100> 

Lr  ar*-  J  W 


=  0 


(38) 


<V100^ 


Introducing  a  time-averaged  stream  function 
ch  sati 

a^ioo^ 


<*'00>  which  satisfies 


1W 


(39) 
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We  may  express  the  above  as 


<V 


u/  >• 

*100 


> 


(40) 


since  any  derivatives  of  <  ^lOO  >  with  regard  to  tj>  are 
zero;  the  biharmonic  equation  is  originally  due  to 
Rayleigh  (ibid.,  4  lines  above  his  equation  (23))  and  de¬ 
rived  anew  by  Riley  (ibid)  by  matched  asymptotic  expan¬ 
sions  , 

The  foregoing  illustrates  the  fact  that  higher- 
order  series  terms  in  expansion  scheme  are  needed,  in 
general,  to  determine  the  present  (1,0,0)  order  stream¬ 
ing  in  the  core;  this  presents  the  entanglement  of  the 
formidable  backward  linking  between  terms  in  series , 
as.  opposed  to  more  direct  forward  linking. 

The  presence  of  the  backward  linking  is  explicitly 
brought  out  here,  since  several  investigators,  without 
recourse  to  the  apparatus  of  the  matched  asymptotic  ex¬ 
pansion,  seem  to  have  unfortunately  erred  on  this  score 


*  To  avoid  any  false  misconception,  we  take  heed  of  the  fact  that 
although  the  streaming  satisfies  the  biharmonic  equation, 
it  is  not  the  Stokes  flow  in  the  sense  that  the  inertia  terms 
are  nil; instead,  the  inertia  terms  are  comprised  of  lower 
order  terms  as  seen  from  equation  (37) . 


by  matching  the  exterior  streaming,  obtained  from  the 

inviscid  equations  expressed  only  to  the  second  order 

cf  disturbances , with  the  one  derived  from  viscous 

consideration  near  the  solid  boundary. 

We  turn  now  to  the  case  or  R  »  1.  The 

s 

fundamental  equation  governing  the  exterior  streaming 
for  this  was  established  by  Riley  (ibid.)  for  incom¬ 
pressible  flow  in  an  otherwise  state  of  rest;  it  cor¬ 
responds  to  the  full  Navier-Stokes  equations  for  steady 

flow,  where  R  replaces  the  Reynolds  number  and  plays 
s 

its  role.  Thus  in  order  to  conquer  it,  nothing 
short  of  wholly  numerical  computation  (Duck  and  Smith, 
1979;  Haddon  and  Riley,  1979)  seem  to  be  effective. 

In  view  of  these  complications,  we  shall 
not  attempt  to  obtain  streaming  in  the  core,  being 
content  instead  that  the  phenomena  of  our  interest 
can  be  explained,  as  seen  shortly,  by  inspecting 
the  behavorial  change  in  streaming  at  the  outer  edge 
of  the  middle  deck,  whose  lead  the  streaming  in  the 
core  follows. 
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5 . 5  Others 

Although  the  remainder  of  the  series  will  not 
contribute  to  the  leading  term  of  acoustic  streaming,  we 
record  here,  for  the  sake  of  completeness,  their  main 
features . 

For  the  time-variant  part  of  (1,0,0)  order,  the 
situation  is  different  from  the  preceeding  streaming  part; 
that  is,  contrary  to  the  latter,  the  former  in  the  core 
is  determined  completely  by  the  (0,0,0)term  in  the  core  and  it  is 
of  inviscid  and  non-heat  conducting  type;  it  is  indepen¬ 
dent  of  its  counterpart  in  the  middle  deck  or  the  time-variant 
part  of  (1,0,0)  order  which  is,  in  turn,  dependent  on  the 
time -variant  part  of  (1,0,0)  order  in  the  core. 

(0,1,0)  series  in  the  core  is  affected  directly 
by  (0,0,0)  order  in  the  middle  deck  and  independent  of 
its  own  counterpart  in  the  middle  deck;  the  boundary 
conditions  represent  an  instantaneous  displacement  effect 
for  the  unsteady  boundary  layer  and  this  turns  out  to  be 
non-zero . 

For  (0,1,0)  order  in  the  middle  deck,  the 
pressure  gradient  in  the  direction  normal  to  the  boundary 
layer  is  no  longer  equal  to  zero,  it  being  balanced  in¬ 
stead  by  a  centrifugal  force  related  to  the  local  curva¬ 
ture;  the  second  coefficient  of  viscosity  is  still  missing 
therein  and  all  these  features  are  analogous  to  the 
second-order  theory  of  the  steady  boundary  layer  (Van  Dyke 
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1962  a,  ibid.)-  In  addition,  the  influence  of  radial 
variation  in  steady  flow  field  appears  for  the  first 
time  in  the  unsteady  boundary  layer  equations. 

All  these,  including  (0,0,1)  series,  possess 
sinusoidal  time -dependence s ;hence ,  once  time -averaged , 
they  all  vanish,  leaving  the  streaming  part  of  (1,0,0) 
order  as  the  leading  term. 
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6.  Solutions  in  the  Middle  Deck 


We  now  present  solutions  in  the  middle  deck. 

Consider  at  first  the  case  corresponding  to  the  thermal 
boundary  condition  where  the  fluctuation  of  temperature  is 
maintained  to  be  zero  on  the  wall. 

From  equation  (34) ,  by  straightforward  calculation 
solutions  for  (0,0,0)  order  in  the  middle  deck  follow 
at  once: 


u000  im  V000 


(r*=l 


■.r  /kr.\  2")  r 

I  6  l1  +  (t)  J  L1  - ex?  (-4n> 


exp  (  -  Pr 


i/2 


,  m  )  J  (  ,  (41-a) 


v000  =  V000  (r*=l>[l  "  exP  (-An)]  •  (41-b) 

^000  =  IT  vooo<riV==1>  [X  "  exP  <"An>]  *  (41_c) 

P000  W***1^1  +  <Y-  1)  exp  (-Pr1/2An)]  ,  (41-d) 

Pooo  -  -  s  <41-c> 

5000  =  "  ri  <Y"1)  V000(r*'1)[1  •  exP  Pr172*")]  •  <41'f) 

6X 


“000  ■  -  [a?l  ,  AH5  <Y"  l)  V000(r**1>[1  -  ex“  <-Pr1/24n)], 

L  J  i=  i  ex 


(41-g) 
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1 

where  A  =  ^ —  (  1  +  i)  and  +  corresponds  to  f^>  0  or 
f^<  0,  respectively. 

Our  principal  formula  ,  the  streamings  at  the 
outer  edge  of  the  middle  deck  near  the  inner  cylinder 
is  immediately  obtained  from  equation  (35-b) .  To 
render  the  final  expressions  explicit,  we  need, in  eIqqq 
terms  ,the  specific  relationship-  between  viscosity  and 
the  temperature.  For  this  we  choose  the  Sutherland's 
formula: 

*  f  ® 

y  =  y  ^0* 

where  y*  is  the  viscosity  at  the  reference  temperature  0*, 
and  s^,is  114°K  for  air.  Then,  the  tangential  streaming 
at  n->«>  or  at  the  outer  edge  of  the  middle  deck  near  the 
inner  cylinder  is  given,  with  the  original  r  restored  in 
place  of  r*,  as 


j 


3/2  e*  + 


"0  +  si 


x 


/ 


3 

1 


2*1 


-  G  + 


(y-D  g 
PtTPr+T)b 


9ex(r’ri> 

slJ 


(«-a) 


where  G 


rl(M  '  Vei£r'rl>~  kWo^r'ri> 


mAex(r-ri) 


This  we  regard  as  one  of  our  central  results . 


(42-b) 


.i 
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Here,  once  again,  we  recall  the  following  meaning 
of  notations  appearing  in  the  above:  u  VQ00 (r=ri)  is 
the  amplitude,  on  the  wall  r  =  r^,  of  the  tangential 
component  of  "inviscid"  linearized  disturbance,  their 
wave  structure  being  in  the  form  of  e^m<^  +  a)t^  ; 

V  ,W  ,A  and  0  are  the  steady  "inviscid"  tangential, 
axial  velocity,  sound  speed,  and  temperature  respectively, 
all  evaluated  on  the  inner  wall;  Pr  is  the  Prandl  number 
Y,  the  ratio  of  specific  heats  and  s^,114°K. 

From  equation  (35-c) ,  the  axial  streaming  is  given 
likewise  by 


[ 


u  V 


000 


<r-ri]  2k 


T  ~  Vex(r=r.)  i 

w_m  - -  _  kWe£=r.)l 

L  i  -* 


j  I  f1  +  (ir)  .  "  G  +  (Y_1)  w$z+i) 


/.  -i  \  pr  n  r 

'  Pr‘  +  1  G  [  2 


<Ur-ri>  H 

«Ur!*ri>  +  =1-1! 

(43) 


The  streaming  in  the  tangential  and  axial  direction 
is  observed  to  be  independent  of  viscosity,  as  to  be 
expected.  To  our  end,  we  need  not  correct  for  the  Stokes 
drift.  We  also  emphasize  that  the  expression  is  not  re¬ 
stricted  to  any  specific,  radial  profile  of  the  steady 
•flow  field. 


.'tf 
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The  radial  streaming  is  given,  from  equation  (35-a), 


<o 

i 


(r=ri)|ai- 


mV 


e 

r 


>£r-ri> 


uV000^r*ri^1  2 
2Aex<r-ri)  J 


X 


+  (y-1) 


where  N  ^r=r^)= 


MejEr-rl> 

RJr'ri> 


(44) 


Nex(r=r^)  being  the  steady  "inviscid"  kinematic  viscosity  on 
the  inner  wall;  the  radial  streaming  is  now  dependent  upon 
viscosity . 

For  the  thermally  insulated  wall,  the  results  can 
formally  be  obtained  from  the  above  by  letting  Pr  ■+*> 

That  is,  for  tangential  streaming,  we  have 

<v*n,  ,  [u  ''ooo(l~ri)1  2m _ 

1  2ri[“-ET  V^-rl)-kwex^i)] 


(Y-1) 


0ex(r-ri)+  S1 


(45) 


and  the  like. 

The  foregoing  expression (42)  or  (45)  for  the  tangen¬ 
tial  streaming  embraces  the  well-known  results  as  its 
special  cases.  In  the  present  representation  of  disturbances, 
ei(m4>+k^  “>t)  *  let  k  =  o  and  in  the  place  of  <j>  ,  introduce 


65 


the  circumferential  distance  x  measured  along  the 
perimeter  of  the  cylinder.  Then  the  disturbance  may 
be  regarded  as  a  plane  wave  traveling  in  the  x  direction  ; 
that  is,  where  l  =  m^ri.  Moreover,  con¬ 

sider  the  fluid  otherwise  in  a  state  of  rest;  then, 

V  =  W  =0  and  A  and  0  are  uniform.  If  the  fluid 
ex  ex  ex  ex 

is  incompressible,  G  vanishes  by  taking  A  -*•  °°  and 
equation  (42)  or  (45)  condenses  to 

<V">  =  I  (u  V000)2  t  .  (46 


This  is  but  a  classical  result  of  streaming  at  the  outer 
edge  of  an  unsteady  shear  layer  (e.g.  Bachelor,  ibid. ,p . 360). 
On  the  other  hand,  if  the  fluid  is  compressible,  G  is 
easily  shown  to  be  unity  and  we  obtain,  for  example, 
from  equation  (45)  for  the  thermally  insulated  wall 


<v'>  = 


V000^ 


20 


1  -  (Y-l)  l 3  - 


ex 


Gex+ 


(47-a) 


where  in  A^and  0^,  the  spatial  dependence  is  now  suppressed. 
The  second  term  in  the  curly  bracket  corresponds  to  the 
streaming  induced  by  the  Reynolds  stresses  which  are 
caused  by  the  variation  of  viscosity  in  response  to  tem¬ 
perature  fluctuation,  coupled  with  the  unsteady  change  in 
strain.  Numerically,  this  is  by  no  means  small;  for 
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instance,  at  0  =  20°C  =  293°K,  it  is  equal  to  0.62  , 
ex 

as  compared  to  unity  of  the  first  term.  This  notwith¬ 
standing,  if  we  choose  to  ignore  it,  the  above  becomes 


<  v'> 


voo(P  2 


Ta 


ex 


(47 -b) 


which  is  also  a  classical  result  for  compressible  flow 
(e.g.  Lighthill , 1978  b,  p.  347). 

Returning  to  the  equation  (42)  or  (45)  for  the  tan¬ 
gential  streaming  in  swirling  flow,  of  utmost  importance 
to  our  objective  is  the  presence  of  the  Doppler  effect  in 
the  form  of 

"  -  ^  Veir-ri>  -  Wex^i> 

in  the  denominantor .  Surely,  then,  this  reveals  that  the 
tangential  streaming  can  reverse  its  direction  around 
turning  points  defined  by 

“  -  rT  V«(r“ri)'  kWex<r'ri>  *  0  ■  <48) 

And  at  the  same  time,  the  absolute  magnitude  of  streaming 

tends  to  beccme  increased  as  such  a  point  is  approached;  the  axial 
streaming  also  exhibits  the  same  features  while  the  radial 
streaming  does  not. 

We  proceed  to  discuss  this  turning  point  in  more 
detail  in  the  next  section. 


7 .  Turning  Points,  Threshold  Swirl 
and  Energy  Separation 


We  define  the  threshold  swirl  as  the  steady  tangential 

velocity  which  corresponds  to  the  turning  point;  that  is, 

V  (r=r.)  which  satisfies  equation  (48).  This  is  deter - 
ex  1 

mined  by  the  competition  between  frequency,  axial  velocity 
and  wave  numbers  of  linear  disturbances.  They  are  inter¬ 
woven  with  each  other  by  the  boundary  condition  for  linear 
waves,  equation  (10-g) .  The  explicit  representations  of 
such  relationships  demand  the  complete  specification  of  the 
radial  profile  of  swirl,  their  derivation  for  a  given 
swirl  being  offered  in  the  Appendix.  By  making  reference 
to  them,  first  we  discuss  the  threshold  swirl  where  the 
radial  steady  profile  is  of  free  vortex  type  confined 
within  an  annulus;  then  we  consider  Rankine  and  forced 
vortices  occupying  a  single  tube. 
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7 . 1  Free  Vortex  between  Co-annular  Cylinders 

Here  our  "inviscid"  velocity  profile  takes 
the  following  form: 

UQ  =  0,  (49-a) 

V0Cr)  -  I  _  (49-b) 

Wq  =  constant,  (49-c) 

for  r.<  r  <r  The  last  condition  satisfies  the  radial 

i  o 

uniformity  of  stagnation  enthalpy,  as  discussed  in 
Section  3.1. 

The  relationship  between  the  frequency  and  the 
above  inviscid  velocity  profile,  an  aforementioned 
prequisite  for  the  determination  of  the  threshold  swirl, 
is  obtainable  numerically  in  the  manner  described  in  (a) 
of  the  Appendix.  Let 


^  'r0)ri' 


''n 


Aex(r*ro) 


(50) 


(51) 


K 


(52) 


The  first  is  the  swirl  Mach  number  referred  to  the  "inviscid" 
acoustic  speed  at  the  outer  cylinder  radius,  A^Cr^r^);  the 
second,  a  dimensionless  frequency  parameter  where  the 
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suffix  n  stands  for  the  order  of  the  radial  mode.  Then, 
illustrated  in  Figure  4  (a)  -  (d)  as  the  solid  lines,  are  the  thus 
computed  curves  of  Xq  corresponding  to  the  lowest  eigen¬ 
value  or  the  first  radial  mode,  drawn  as  functions  of  M, 
for  the  various  radius  ratios  X  ;  the  tangential  modes  corres¬ 
pond  to  the  first  (m=l)  and  the  second  (m=2)  andK  =  0,l;the  ones 
lying  in  the  range  of  positive  correspond  to  the  dis¬ 

turbances  spinning  in  the  same  direction  as  the  steady 
swirl;  those  for  negative  correspond  to  the  distur¬ 

bances  spinning  in  the  opposite  direction.*  Observe  that 
over  a  wide  range  the  frequency  is  nearly  proportional  to 
the  swirl  Mach  number,  a  feature  pointed  out  in  Section 
1.1  in  connection  with  the  annular  cascade  data.  The 
virtual  linearity  suggests  immediately  an  analytical  so¬ 
lution  by  resorting  to  the  expansion  in  terms  of  the  swirl 
Mach  number;  such  an  approximate  formula  is  given  by 
equation  (A-8)  of  the  Appendix,  shown  as  chain  lines  in 
Figure  4,  and,  on  the  whole,  the  results  stand  in  reasonable 
agreement  with  the  wholly  numerical  results  .  The  agreement 
is  closer  for  a  narrow  annulus;  this  will  afford  us  to  use 
and  exploit  the  analytical  results  for  such  a  case,  as  to 
be  seen  shortly. 

With  the  aid  of  these  relationships,  equation 
(42)  or  (45)  furnish  the  result  for  tangential  streaming. 
Figure  5  (a)  -(d)  typify  the  one  for  several  radius  ratios;  in  all 
of  them,  the  dimensionless  tangential  streaming  near  the  inner  cylinder 

One  can  also  represent  this  backward  traveling  wave  by  the 
one  with  negative  in  in  the  positive  range  of  M^  . 


is  plotted  as  a  function  of  the  steady  swirl  Mach  number,  ; 
the  mode  in  the  radial  direction  is  the  fundamental, 
m  “  1,  and  k=  0,  and  the  frequency  relationship  used 
therein  is  that  obtained  by  the  numerical  procedure; 
the  solid  lines  correspond  to  the  thermal  condition 
where  the  fluctuation  of  the  temperature  is  main¬ 
tained  to  be  zero  on  the  wall  of  the  inner  cylinder, 
the  broken  lines  to  thermally  insulated  walls.  There, 
the  positive  value  of  streaming  means  that  its  direction 
is  the  same  as  that  of  steady  swirl,  negative  being 
obviously  the  other  way.  Thus,  the  figures  indeed 
display  the  reversal  of  the  streaming  direction  at 
the  threshold  swirl  Mach  number.  And  the  absolute 
magnitude  of  streaming  tends  to  become  sizable  in  its 
neighborhood  with  the  threshold  swirls  serving  as 
asymptotes.  To  obtain  the  total  magnitude  of  "d.c." 
swirl,  we  only  add  or  subtract  the  streaming  to  the 
steady  swirl  imposed,  of  course;  therefore,  beyond  the 
threshold,  the  entire  swirl  becomes  reduced  near  the 
inner  cylinder.  Notice  also  that  the  threshold  swirl 
Mach  number  decreases  as  the  ratio  of  the  outer  to 
inner  radius  of  the  cylinder  \  is  progressively  increased. 

The  point  appears  more  directly  in  Figure  6  (a) -(b) 
where  the  threshold  swirl  Mach  numer  of  the  foregoing 
first  radial  mode  is  drawn  against  the  radius  ratio 
for  the  first  and  second  tangential  modes,  m  -  1  and 
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m  “  2;  as  an  example,  for  the  radius  ratio  of  5,  the 
threshold  swirl  for  m  ■  1  and  ic  ■  0  is  only  0.37, 
a  low  value  for  tangential  velocity  indeed. 

The  threshold  swirl  is  independent  of  the 
thermal  condition  at  the  wall;  it  corresponds  to  the 
disturbances  spinning  in  the  same  direction  as  the 
steady  swirl;  no  threshold  is  found  for  the  ones  in 
the  opposite  direction.  In  addition,  we  can  show  that 
for  any  values  of  k,  whether  negative  or  positive,  the 
lowest  threshold  swirl  Mach  number  occurs  at  k  =  0. 

Figure  6  also  reveals  that  as  the  radius  ratio  X 
approaches  unity,  the  threshold  swirl  becomes  exceed¬ 
ingly  high.  In  fact  we  can  prove  that  for  such  a  narrow 
annulus,  in  its  limit,  the  threshold  swirl  approaches 
infinity;  that  is,  whatever  the  amount  of  swirl,  no 
retrogression  of  the  tangential  streaming  occurs. 

This  comes  from  the  fact  that  as  exhibited  by  the  afore¬ 
mentioned  analytical  solution  (A  -  19)  of  the  Appendix, 
the  frequency -swirl  relation  for  the  lowest  eigenvalue 
and  tc=  0  becomes  such  that  it  corresponds,  in  effect, 
to  the  plane  wave  --  traveling  with  the  phase  speed 
equal  to  the  sum  and  difference  of  swirl  and  ^(r  “ro^  * 

To  the  extent  no  reversal  of  the  streaming  direction  takes 
place  for  the  plane  wave,  no  retrogression  of  tangential 
swirl  occurs  in  the  present  limit.  As  a  matter  of  fact, 
we  can  expressly  prove  that  for  such  a  case,  the  expression 
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of  equation  (42)  or  (45)  is  reduced  to  such  as  equation 
(47);  the  denominator  now  becomes  A^,  a  constant,  and 
therefore  no  reversal  takes  effect . 

For  the  foregoing  reason,  the  reversal  in  the 
streaming  direction  is  intimately  connected  to  the  pre¬ 
sent  three-dimensional  geometry,  which  compels  the 
lowest  mode  of  disturbances  to  deviate  from  the  plane 
waves ;  also,  no  less  indispensable  is  the  presence  of 
steady  swirl. 

For  the  outer  cylinder,  Figure  7  shows  that 
the  tangential  streaming  at  the  outer  edge  of  its  middle 
deck.  Although,  contrary  to  Figure  5,  no  reversal  is 
experienced  therein,  this  is  still  indirectly  affected 
by  the  significantchange  of  streaming  near  the  inner 
cylinder:  the  reduction  in  tangential  streaming  of  the  latter 
induces,  in  the  annular  core,  a  complex  radial  redistri¬ 
bution  of  the  angular  momentum,  whose  precise  details  are 
outside  the  scope  of  the  present  work;  however,  d\  e  to  the 
fact  that  the  torque  acting  upon  the  annular  core  remains  vir¬ 
tually  zero,  we  can  say  that  this  decrease  has  to  be  com¬ 
pensated  by  the  corresponding  increase  of  tangential  velocity 
elsewhere, including  the  region  near  the  outer  cylinder. 

When  synthesized,  the  transformation  of  the  total 
d.c.  swirl  above  the  threshold  appears  as  shown  schematically 
in  Figure  8;  thus,  post-conversion  behavior  would  plainly 
separate  the  total  temperature  in  the  radial  direction, 
even  if  the  static  temperature  remained  uniform  (we  shall 


discuss  this  latter  aspect  in  the  next  subsection  7.2  . 

The  almost  linear  dependence  of  frequency  upon 
swirl,  the  existence  of  threshold  swirl,  the  reduction 
of  total  d.c.  component  of  swirl  near  the  inner  cylin¬ 
der  beyond  the  threshold  and  the  radial  gradient  of 

total  temperature  or  the  Ranque-Hilsch  effect -  all 

appear  to  agree,  at  least  qualitatively,  with  the  key 
features  mentioned  on  the  annular  cascade  in  Section  1,1. 
(For  now,  further  quantitative  comparison  is  beyond  us 
for  want  of  more  data  surrounding  the  threshold  swirl.) 

In  reference  to  Figure  5,  as  the  threshold 
swirl  is  approached  with  increasing  amounts  of  streaming, 
the  assumption  of  small  parameter  of  a  and  6  breaks  down; 
moreover,  the  change  induced  in  the  base  flow  pattern 
starts  to  modify  the  frequency-swirlrelationship.  Neverthe¬ 
less  ,  we  assert  that  the  phenomenon  is ,  at  the  very 
least,  an  indication  of  what  is  to  be  expected  of  the  actual 
one  at  or  near  such  a  condition,  the  only  substantive 
modification  in  reality  being  more  gradual  reversal  of 
the  streaming  direction  rather  than  the  abrupt  one  across 
the  asymptotes  depicted  in  the  figure. 

Still  referring  to  Figure  5,  just  below  the 
threshold  swirl,  we  observe  that  the  streaming  enlarges 
its  magnitude;  this  would  increase  the  total  d.c.  component 
of  swirl  near  the  inner  cylinder.  We  venture  to  state, 
however,  that  in  reality  such  development  is  unlikely  to 
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occur  due  to  the  ensuing  unbalance  of  radial  pressure 
gradient  and  centrifugal  force.  This  expels  outward  the 
lumps  of  fluid,  which  is,  at  the  same  time,  swept  down¬ 
stream.  It  is  replaced  afresh  by  another  fluid 
element  replenished  from  the  mainstream.  This  new  lump, 
while  migrating  inward  through  the  virtually  inviscid 
region,  increases  its  tangential  velocity  in  accord  with 
Kelvin's  theorem  until  it  attains,  upon  arrival  at  the 
inner  radius,  the  originally  imposed  steady  tangential 

i 

velocity;  thereabout  the  steady  swirl  remains,  in  effect, 
unchanged . 

On  the  other  hand,  above  the  threshold  swirl, 
the  reduction  in  total  d.c.  component  of  the  swirl 
can  occur  indeed.  The  unbalance  of  radial  pressure 
gradient  and  centrifugal  force  tends,  in  this  instance, 
to  submerge  the  lump  of  fluid  toward  the  wall  of  the 
inner  cylinder.  But,  it  is  prevented  from  doing  so  by  the 
radial  outflow  resulting  from  the  steady  radial  velocity 
ejected  from  the  steady  boundary  layer, in  addition  to  the 
radial  component  of  acoustic  streaming  directed  also  and 
always  outward,  as  seen  from  equation  (44);  hence,  the 
fluid  element  is  forced  to  remain  adhered  to  the  inner 
cylinder  with  its  magnitude  of  tangential  velocity  now 
reduced . 
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7.2  Rankine  Vortex 

We  cons.ider  the  Rankine  vortex  within  a 
single  tube  represented  by 

VQ(r)  =  flr  ,  for  0  <  r  <  r±,  (53-a) 

VQ(r)  “  r  »  for  r±  <  r  <  rQ,  (53-b) 


2 

where  r  =  fir^  and  r^  now  denotes  the  radius  at  tne  interface 
of  a  forced  and  free  vortex.  In  addition,  UQ  =  0  and 
W^(r)  is  an  arbitrary  function  of  r. 

Furthermore ,  we  are  interested  in  the  unsteady 
disturbance  with  k  »  0  or  ic  »  0  and  the  first  radial  mode. 

As  shown  in  equation  (A- 24  b)  of  the  Appendix,  the 
frequency -swirl  relationship  for  such  a  case  can  adequately 
be  represented  by 


-2|m| 

u  =K|m|  -  1  +  X 


(54) 


+  sign  corresponding  to  m  >  0  or  m  <  0  respectively, 
provided  the  swirl  Mach  number  at  the  interface  remains 
well  in  the  subsonic  range.  This  formula  asserts  that 

« 

the  frequency  is  proportional  to  the  swirl,  a  distinguishing 
characteristic  of  vortex  whistle  first  discovered  by 
Vonnegut  (1954);  see  also  Chanaud  (1963,  1965).  Its 
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linearity  is  to  such  an  extent  and  with  such  accuracy 
that,  by  creating  swirl  and  measuring  its  frequency, 
it  can  be  exploited  for  use  either  as  a  sensor  of 
aircraft  speed  (Nichlas,  1957)  or  a  flow  meter,  provided 
in  the  latter  the  ratio  of  swirl  to  axial  velocity  is 
made  to  remain  constant  (Rodely,  Chanaud  and  White, 

1965) .  (Equation  (54)  is  applicable  even  to  incompres¬ 
sible  flow,  as  may  be  inferred  to  be  so  by  noting  the 
absence  of  acoustic  speed  and  density  thereof;  this  can 
also  be  directly  verified  by  examining  the  incompressible 
solution  of  Kelvin  (ibid . ) .  The  proportionality  of  surge 
frequency  to  rotational  speed  of  water  turbines,  as 
observed  earlier  by  Rheingaus  (1940)  in  draft  tubes  of 
hydroelectric  plants  is  yet  another  practical  manifestation 
of  the  vortex  whistle  in  incompressible  flow.) 

Although  both  for  a  Rankine  and  free 
vortex,  the  relationship  between  the  frequency  and  swirl 
is  found  to  be  linear,  the  difference  between  full 
proportionality  (i.e.  id  -*■  0  as  0  -*■  0)  for  the  present 
Rankine  vortex  and  mere  linearity  (id  /  0  as  r  -*■  0) 
for  the  free  vortex  of  the  preceeding  section  is  enough 
to  compel  a  notable  change  in  the  turning  point,  as  will 
be  seen  below. 

Making  use  of  (54) ,  the  tangential  streaming 
can  at  once  be  obtained  from  either  equation  (42)  or  (45) 


by  replacing  the  subscript  i  with  o.  Figure  9 
exemplifies  such  calculated  results;  the  tangential  streaming 
at  the  tube  periphery  is  rendered  into  dimensionless  form  and  drawn 
as  a  continuous  function  of  circumferential  wave  number 
m,  though,  in’  reality,  m  of  course  takes  only  integer 
values;  the  swirl  Mach  number  at  the  interface 
r/r.A  (r=r  ) , is  taken  to  be  0.2,  X  being  2.  Observe 
that  the  direction  of  the  tangential  streaming  is 
always  positive  or  in  the  same  direction  as  the  steady 
swirl;  this  is  so, not  only  for  positive  values  of  m 
corresponding  to  the  disturbances  spinning  in  the  same 
direction  as  steady  swirl,  but  also  for  negative 
m  corresponding  to  the  ones  in  the  opposite  direction. 
Observe  also  that  the  magnitude  of  streaming  increases 
as  the  circumferential  wave  number  is  decreased.  In 
fact,  at  m  =  +  1  or  the  first  tangential  mode,  it  can  be 
readily  shown  from  equations  (42)  or  (45)  combined  with 
(54)  that  the  magnitude  of  steady  streaming  becomes  in¬ 
finitely  large,  regardless  of  the  value  for  steady  swirl, 
or  of  Xor  of  the  thermal  condition  on  the  wall;  this  implies 
that  the  turning  point,  contrary  to  the  free  vortex,  does 
not  exist  for  the  present  case.  (Of  course,  in  reality, 
the  magnitude  does  not  become  that  much,  but  again  we 
take  it  as  an  indication  of  what  actually  happens.) 

For  the  reason  stated  in  the  foregoing  sub¬ 
section  7.1  ,  the  unbalance  between  the  radial  pressure 


gradient  and  the  centrifugal  force  work  in  such  a 
way  as  to  bring  the  increase  of  streaming  at  the  tube 
perimeter  into  being. 

Schematically  drawn  in  Figure  10  is  the 
conversion  of  the  original  Rankine  vortex  to  a  forced 
one  by  streaming.  Acoustic  streaming,  if  induced, 
being  present  regardless  of  the  axial  position  in  the 
tube,  the  transformation  can  immediately  take  effect 
right  at  its  entrance.  This  metamorphosis  toward 
a  forced  vortex  and  this  alone,  to  say  nothing  about 
the  static  temperature  gradient,  again  tends  to  separate 
the  total  temperature  into  a  colder  stream  near  the  center  and 
the  hotter  stream  near  the  perphery  of  the  cylinder. 

If  and  when  the  geometry  of  the  tube  arrangement 
favors  the  excitation  of  such  particular  unsteady  dis¬ 
turbance  as  the  one  with  the  lowest  eigenvalues  in  the 
radial  direction,  m  =  1,  and  k  =  0,  a  forced  vortex 
filling  out  the  entire  tube  is  always  formed  and  the 
total  temperature  becomes  separated,  even  near  the  inlet; 
otherwise,  the  Rankine  vortex  and  the  total  temperature 
remain  unchanged.  This  appears  to  explain  the  dichotomy 
of  the  redial  profile  mentioned  in  the  Introduction:  a 
forced  vortex  for  a  Ranque-Hilsch  tube  with  separated 
total  temperature  and  Rankine  vortex  for  others  without 
any  such  separation. 
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By  appealing  to  the  streaming  as  the  predominant 
mechanism,  we  can  also  readily  explain  the  other  character¬ 
istics  of  the  Ranque-Hilsch  effect:  the  radial  difference 
of  total  temperature  is  known  to  increase  as  the  pressure  of 
the  incoming  air  is  raised.  As  the  pressure  is  stepped  up,  the 
amplitude  of  disturbances  as  well  as  the  amount  of  steady 
swirl  increases  (Cassidy  and  Falvey,  1970);  from  equation 
(42-a) ,  this, in  turn  , is  observed  to  induce  more  streaming 
and  hence  more  energy  separations . 

So  far,  we  have  not  touched  upon  the  role  of 
static  temperature.  The  steady  state  or  time-averaged  component 
of  static  tenperature  induced  by  the  unsteady  disturbances  at  the  outer 
edge  of  the  middle  deck  is  always  of  insignificant  magni¬ 
tude;  hence  it  would  not  sensibly  affect  the  basic  radial 
"inviscld"  pattern  of  static  temperature,  determined  uniquely 
by  swirl,  for  any  pattern  of  which  the  following  holds: 
higher  temperature  at  the  outer  radius  and  a  colder  one 
near  the  center.  In  fact,  the  data  of  Eckert  and  Hartnett 
and  Lay  (ibid.)  all  show  such  radial  distribution  of  static 
temperature  for  Ranque-Hilsch  tubes;  the  actual  preservation  of  this  in 
the  presence  of  unsteadiness  existing  in  the  flow  may 
indeed  involve  the  Knoernschild  effect  (ibid.),  which 
is  applicable  whether  the  disturbance  is  of  organized  origin 
or  of  stochastic  nature.  In  any  event,  this  difference  in 
static  temperature,  when  added,  aids  to  separate  furthermore 
the  temperature. 


2* 
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Once  the  formation  of  a  forced  vortex  is 
thus  predicted,  then,  based  on  this,  together  with 
assumed  distribution  of  static  temperature  such  as  an 
adiabatic  one,  and  from  these  two  alone,  it  is  a  simple 
matter  to  construct  the  so-called  performance  curve 
resembling  the  one  given  by  Hilsch;  in  other  words,  upon 
only  the  above  two  the  performance  curve  rests;  hence, 
once  they  are  somehow  gotten  in  any  theory  whatsoever, 
the  performance  curve  follows  at  once  and,  contrary  to 
what  might  be  believed,  the  mere  reproduction  od  the 
curve  does  not  substantiate  the  postulated  mechanism. 


.  ■  -!X  .  - 
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7.3.  Forced  Vortex 

Consider  a  forced  vortex  given  by 

Vq  (r)  =  flr  ,  0  <  r  <  rQ  ,  (55) 

and  Uq  =  Wq  =  0.  With  the  aid  of  the  relationship 
between  w  and  ft  obtained  numerically  by  Sozou  (1969) , 
we  can  show  that  no  turning  point  exists  for  such  swirl. 


In  the  foregoing,  we  have  considered  situations 
where,  in  addition  to  the  steady  swirl  and  axial  flow, 
both  the  tangential  and  axial  wave  numbers  are  specified 
and  the  frequency  is  to  be  determined;  on  the  other  hand, 
if  the  frequency  of  the  disturbances  is  enforced  upon 
the  swirling  flow,  held  fixed  and  the  axial  wave  num¬ 
ber  is, to  be  determined,  such  a  case  may  be  treated 
in  a  similar  manner,  the  streaming  still  being  given  by 
equations  (42)  to  (45) . 
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8 .  Summary 

t 

In  this  last  section,  we  summarize  the  foregoing 
results  on  streaming,  interpreting  them  now  in  physical 
terms  from  the  outset. 

We  begin  by  recalling  that  for  a  wave  traveling  uni- 
directionally  through  a  gas  in  an  otherwise  state  of 
rest,  the  effect  of  viscosity  manifests  itself  in  the  form 
of  the  tilted  trajectory  of  a  fluid  particle  near  the 
solid  wall;  it  is  this  skewness  that  gives  rise  to  the 
streaming.  By  skewness  is  meant  that  the  axes  of  the 
particle  path  tracing  a  closed  loop  are  slanted  as  re¬ 
ferred,  for  example,  to  the  normal  to  the  wall  surface. 

We  draw  the  control  surface,  parallel  to  the  wall  and 
intersecting  the  loop,  and  reckon  the  momentum  trans¬ 
ported  by  the  lump  of  fluid  circling  along  the  loop,  in 
and  out  of  the  surface.  Then,  the  very  asymmetry  of 
slanting  trajectory  results  in  the  unbalance  of  the  mo¬ 
mentum  budget,  once  averaged  temporally;  the  Reynolds 
stresses,  thus  generated  and  of  orderly  origin,  act  upon 
the  fluid  layer  to  set  forth  the  streaming. 

If  the  wave  is  propagating  in  the  x  direction  with 
waveform  given  as  ue^*1*  “  at  the  outer  edge  of  the 
unsteady  boundary  layer,  the  streaming  in  the  x  direction 
at  the  edge  takes  the  form  of 

~2 

S-i  (  ) 

when  by  positive  value  we  mean  that  its  direction  is  the 
same  as  that  for  the  wave  propagation. 
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If  the  wave  is  spinning  in  the  circumferential 
direction  with  waveform,  given  by  ats 

the  outer  edge  of  the  unsteady  boundary  layer  formed 
over  a  cylindrical  surface  of  radius  r^,  the  streaming 
in  the  tangential  direction  follows  immediately  from 
the  above  as 


“2 
u  m 


(...)  . 


If,  in  addition,  the  gas  itself  swirls, with  its 
magnitude  at  the  outer  edge  given  as  Vex,  the  Doppler 
shift  changes  the  above 


so  that  to  an  observer  situated  on  the  frame  of  reference 

rotating  with  peripheral  speed  Vex>  this  would  apparently 

become  the  same  as  the  one  for  the  gas  which  is  otherwise 

at  rest.  The  auove  can  readily  be  generalized  to  the 

three-dimensional  situation  where  the  waveform  is  now 

represented  by  e^m<^  +  wt  ^ ;  besides  swirl,  the  steady 

velocity  now  possesses  an  additional  axial  component  Wgx. 

Then,  in  place  of  the  foregoing  Doppler  shift  term 
V  V 

nv  gy 

w  -  m  — —  ,  we  only  have  to  replace  with  o»-  m  ~  -  kVJ  , 

*  i  i 

the  frame  of  reference  being  now  in  screw  motion,  advancing 


84 


axially  with  W  ,  while  rotating  with  V  .  The  elaborate 
apparatus  of  the  matched  asymtotic  expansion  applied  to  the 
full  Navier-Stokes  equation,  has  led  us,  in  fact,  to  this 
form  as  the  leading  term  of  streaming  in  the  tangential 
direction,  the  entire  formula  of  which  is  given  as  equation  (42)  or  (45) 
The  very  form  of  the  above  Doppler  shift  implies  that 
the  tangential  streaming  can  change  its  rotational  direction. 

V 

For  w  -  m  — —  -  kVJ  >0,  it  rotates  in  the  same  direction 

ri  ex  .  V 

as  the  steady  swirl;  for  oj-  m  -  kW  <0,  the 

r^  ex 

streaming  is  now  rotating  in  the  opposite  direction  - 

retrogression  occurs.  As  such  a  turning  point  is  approach¬ 
ed,  the  absolute  magnitude  of  streaming  increases  sizably. 

The  turning  point  is  determined  by  the  relative  values 
of  competing  factors  such  as  frequency  w ,  steady  flow  field, 

Vex  and  Wex,  and  the  wave  numbers,  m  and  k;  in  turn,  they 
are  dependent  upon  the  specific  radial  profile  of  swirl 
and  the  geometry  of  cylindrical  configurations . 

For  a  free  vortex  distribution  contained  between  two 
concentric  cylinders  of  radii,  r^  and  r^ ,  the  threshold 
swirl  corresponding  to  the  turning  point  decreases  as  the 

ro 

radius  ratio  X  =  —  is  increased,  in  a  fashion  depicted  in 
Figure  6.  Above  the  threshold,  the  streaming,  when  added 
to  the  steady  swirl,  tends  to  reduce  the  total  d.c. 
amount  of  swirl  near  the  inner  cylinder. 
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For  a  Rankine  vortex,  within  a  single  cylinder 
and  subject  to  the  particular  wave  corresponding  to 
the  lowest  eigenvalue,  m  =  1,  and  tc  -  0,  there  is 
no  threshold  swirl;  instead,  any  amount  of  swirl  in¬ 
troduces,  along  the  cylinder  periphery,  a  streaming 
rotating  in  the  same  direction  as  the  steady  swirl 
and  of  considerable  strength,  converting  it  into  a 
forced  vortex.  For  waves  other  than  this,  the  effect 
of  streaming  remains  insignificant. 

Be  it  a  free  vortex  or  a  Rankine  vortex,  if  the 
original  steady  swirl  metamorphoses  into  a  forced 
vortex,  then  this  and  this  alone,  suffices  to  give  rise 
to  the  Ranque-Hilsch  effect  (though the  presence  of  depressed 
static  temperature  distribution  near  the  center  also 
furthers  this) . 

The  foregoing  appears  to  account  for  the  phenomena 
described  in  connection  with  the  annular  cascade  and  the 
issues  raised  for  the  Ranque-Hilsch  tube  alike, as  de¬ 
tailed  in  the  Introduction. 

Finally,  we  close  this  paper  with  the  remark  that 
to  the  extent  the  temperature  separation  arises  due  to 
the  unsteadiness  in  flow,  the  present  subject  is  the 
converse  to  the  phenomena  of  Rijke  tube  and  thermally 
driven  acoustic  oscillation  of  liquid  helium  (e.g.  Clement 
and  Gaffney,  1960)  where  the  very  difference  in  tempera¬ 
ture  gives  rise  to  unsteady  disturbance;  and  in  broader 


86 


contexts,  to  the  extent  that  we  seek  the  present 
mechanism  as  of  organized  origin  distinct  from  the 
stochastic  process,  this  falls  under  the  same  mor¬ 
phological  class  as  the  study  of  large-scale  structure 
in  the  mixing  layer  (Brown  and  Roshko,  1974). 

The  author  is  indebted  to  Mr.  C.  E.  Danforth  for 
calling  his  attention  to  the  phenomenon  of  the  annular 
cascade  and  its  connection  to  the  Ranque-Hilsch  effect; 
to  Dr.  M.  E.  Goldstein,  for  helping  him  to  clarify 
several  points  by  raising  pertinent  questions  ;  to  Drs. 
Caruthers  and  Maus  for  providing  valuable  comments  on 
the  original  manuscript;  to  Drs.  Cassidy,  Falvey,  Flandro, 
Keyes,  Lay,  Ragsdale,  Savino,  Sforzini  and  Vonnegut,  for 
answering  to  his  various  queries  related  to  their  works . 
The  figures  are  based  on  the  computation  performed  by 
Mr.  J.  Q.  Chu  and  Ms.  D.  E.  Gonzalez;  Lt .  J.  M.  McGee 
also  contributed  to  this .  The  work  is  supported  by  the 
Air  Force  Office  of  Scientific  Research  under  Contract  No. 
F  49620-78-C-0045 . 
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This  appendix  concerns  the  relationship  between  the 
frequency  and  the  steady  swirl  corresponding  to  section 
7.1  and  7.2;  it  will  be  obtained  from  the  unsteady,  inviscid 
linearized  approximation,  equation  (10)  or  equivalently 
equation  (33) . 

a)  Free  Vortex  Between  Co-Annular  Cylinders 
For  this,  equation  (49)  specifies  the  steady  profile. 
Choosing  r*v  or  the  moment  of  angular  momentum  as  our  primary 
dependent  variable,  we  can  extract,  from  equation  (10)  and  with  the 
aid  of  equation  (8),  the  following, single, second-order  differential  equation: 


(A-l) 


with  the  boundary  condition  given  by 

-  0  at  r*=  1  and  i*=  X.  (A-2) 

The  definitions  of  ,  xn  and  k  are  given  in  equation  (50) 
to  (52) .  The  eigenvalues  xn  are  to  be  determined  once  the 
values  of  M^,  k  ,  m  and  X  are  given; 
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Figure  4  in  the  text  has  been  obtained  by  the  standard 
eigenvalue  search  procedure  using  the  fourth  order 
Runge-Kutta  method. 

In  addition  to  this  wholly  numerical  procedure. we 
can  also  obtain  an  approximate  analytical  representation 


of  eigenvalues  and  this  we  discuss  next. 


(0) 


Consider  at  first  the  case  of  no  swirl.  Then  r*v  v J 
be  expressed  as  a  linear  combination  of  Bessel  functions 
such  as 


r*v  ^  =  A 


[Jm<zn 


"  Ym(zn 


J'  (O 


Y  U 
m  n 


n'  ' 
)  - 


(A-3) 


where  A  is  a  constant,  zr  is  given  by 

V<zn*>  V  <zn>  -  V<zn*>  V<zn>  "  °- 


where  z  *  is  related  to  z  by 
n  n  ' 


(A-4) 


z  *  =  z  X 
n  n 


(A-5) 


In  (A-4) ,  the  primes  denote  the  derivatives  with  respect 
to  the  argument  of  the  Bessel  function.  zfi  is  related  to 
Xn  by 


2  .  2 
k  +  zn 


(A-6) 
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Based  upon  this,  when  the  swirl  is  present,  we  expand 

* _ 

rv  and  *-n  terms  of  the  swirl  Mach  number,  M.  and  retain 
n  <P 

only  the  linear  term: 


r*v 
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rv 
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(A-8) 


where  r*v^and  Snm  are  t6  be  determined.  Substituting  into 
(A-l) ,  retaining  only  terms  up  to  0  (M^,)  and  satisfying  (A-2) 
and  making  use  of  the  integral  formulae  involving  products  of 
the  Bessel  function  (e.g.  Watson  1966, pp. 134-137) ,  we  obtain, 
for  m  >  0 


(A-9) 


where 


A1  =  2V<zn>V<zn*>  hi™  ~  V(zn>  V<zn*>  *12  <J  > 
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^11’  *12*  etc*  are  S*-ven  by 


IU(JY) 


~  2k  Jo(z)Yo<z)  +  2  5Z  J 

L  £=1 

* 

] 


l<2>  \  (z) 


z  =  z_ 


+  Vz>  V2> 


n 


(a-: 


z  =  z 


n 


I12(j2) 


'  K  [Jo 


m-1 


(z)  2  +  2  £  J*  (z)z  +  J 

«.=! 


m(z>  2J 


z  = 


(A- 


rl3<Y  >  ' 


m-1 


2m  I  Vz>  +  2  2 


1 


\  U)  +  Ym(z) 


■] 


z  = 


z  = 


(A- 


*21  (JY) 


1  ,2  [ 
%  L 


2Jm<z)  VZ>  "  J™-l<z)Y»n4-l(z)  -Jm+l(z)Y: 


(A- 


J22  (J  ) 


|l  z2  [  Jm(z)  2  "  Jm-l(z)  Jm+l(z)  ]jz  «  z*  . 


(A- 


M 


,2S  _  \1  2 !  „  ,_v  2 


j23  a")  =  Ym(z)  -  YJn_1  (z)  Y^Cz) 


](■■■" 
J  ’  z  =  z, 


n 


(A 


From  (A-12)  to  (A-17),  (z)  ]  zIzn  denotes  the  following 
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[fOOj  l~_Zzn*  =  f  <zn*>  "  f<2n>  *  (A"18) 

(If  m<0,  S  remains  the  same  for  m>  0. ) 
nm 

The  values  of  the  lowest  eigenvalue  Zq  and  corresponding  S^, 
computed  from  equation  (A-4)  and  (A-9),are  listed  in 
Table  1. 


X 

_ ^0_ 

m=l 

m=2 

S01 

CM 

o 

M 

1.01 

0.99503 

1.99006 

0.99009 

0.99016 

1.05 

0.97571 

1.95141 

0.95202 

0.95202 

1.1 

0.95274 

1.90550 

0.90772 

0.90770 

1.5 

0.80509 

1.60810 

0.64760 

0.64420 

2.0 

0.67734 

1.34060 

0.45551 

0.43714 

3.0 

0.51362 

0.97749 

0.25366 

0.21202 

5.0 

0.34102 

0.60694 

0.10241 

0.-07354 

7.0 

0.25240 

0.43556 

0.05240 

0.03686 

10.0 

0.18035 

0.30529 

0.02516 

0.01793 

TABLE  1 . 

The  comparison  of  the  foregoing  analytical  repre¬ 
sentation  with  the  numerical  procedure  over  a  wide  range 
of  swirl  is  given  in  Figure  4. 

As  a  special  case  of  the  above,  we  focus  expressly 
our  attention  to  the  narrow  annulus;  let  X-*!. 
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Then,  we  can  show  explicitly  that 


and 


S 


Om 


For  k  =  0,  this  leads  immediately  to 


a)  =  A  (r  =r  ) 
ex'  oJ 


lm_| 

r. 


mr _ 

rori 


from  equations  (A-8)  and  (51).  Then  the  disturbance  repre¬ 
sented  by  e^111^  becomes, for  the  narrow  annulus 


(A-19) 

The  term  within  the  second  £  J  shows  t*hat  this  corresponds  to  a  plane 
wave  propagating  at  the  acoustic  speed  while  being  convected 
at  the  swirl  velocity,  as  to  be  expected. 

If  we  change  our  viewpoint  slightly  and  regard  u  as 
given  and  k  to  be  sought,  then  by  insisting  that  k  is  real  in 
(A-8)  and  (51),  the  so-called  cut-off  condition  for  the 
propagation  of  disturbances  through  ducts  can  be  obtained; 
this  generalizes  the  Tyler-Sofrin  formula  (1962)  by  including 
the  presence  of  both  the  swirl  and  the  axial  flow. 


II  Hill  l  II  HIM  I  ir 


»-  V 
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b)  Rankine  Vortex 

For  the  radial  profile  given  by  equation  (53) ,  the 
boundary  conditions  for  the  unsteady  disturbances  are: 


u  (r  =rQ)  =  0, 


(A-20a) 


p  and  u  are  continuous  at  r  =  r~  (A-20b) 

Sozou  and  Swithenbank  (1969)  calculated  the  corre¬ 
sponding  frequency  relationship  by  wholly  numerical 
computation.  (Their  treatment  is  limited  to  the  two- 
dimensional  flow.  Although,  as  stated  in  Section  7.2,  the 
axial  velocity  is, at  present .allowed  to  exist,  as  long 
as  k  is  set  equal  to  zero,  its  presence  will  affect  only 
the  axial  disturbance,  leaving  the  rest  same  as  two- 
dimensional  flow.)  Here  we  offer  an  approximate  analytical 
representation,  choosing  fT  as  a  primary  dependent  variable. 

In  the  innermost  core  of  a  forced  vortex  region  defined  by 
0<r<r^,the  governing  equation  is  given,  upon  combining 
equation  (10)  together  ,by 


d2  p  +  [~1  +  (2y  -  3)  02  r  1  dp 

Lr  A02(r)  J  ^ 

+  r^a-23.-  +  tty- s.2  .  4  +  »  i  _ ,  n 

LxiAo2(r)  r  Ao  <r>J  p 


(A-21a) 
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where  AQ<r),x^,  and  D  are  given  by  the  following: 

A02<r>  =  Aex2(r"ro)+  Z<Y  ‘  15  “2  <r2-  ro2>- 
=  flm  -  w, 

D  =  (Qm  -to)  2  -  4fi2  . 


(A-21b) 

(A-21c) 

(A-21d) 


u  is  related  to  p  by 

WD  I  'l  ^  +  ? 

.2 


R0(r)D 


2mA0  (r)ft 


-  (2  -y  )  Xx  r  n' 


(A- 21  e) 


On  the  other  hand,  in  the  free  vortex  region  defined 
by  ri<r  <  ro«  the  governing  equation  is  given  by 


,2- 


dr 


Vr> 


Ao 

'  -  3)  r2  4m r 

1 - n -  +  o 

(r)r J  X2(r)rJ 

£- 

'x22 (r) 

_  m2  +  2m(2y-  3) r  3 

8(mr)2 

Ao2(r) 

7  X2(r)r\Z(r) 

X22(r)r6 

4m  £ 

X2(r)r4 

.  2(2-y)  r2 

AQZ(r)r4 

p  = 

0  , 

(A-22a) 

"  Aex2<r 

=r0)  +  ^(T-l>[(!2 

<ri2-r, 

02>t2( 

mr 

W  , 

(A-22b) 

7  ' 

(A-22c) 

nd  u  is  related  to  p  by 
i  )  “0 


An2(r) 


u 


RQ(r)  |  X2 (r) 


d£ 


dr 


+ 


2mA0Z(r)r 


TT 

L  X2  r 


*7^ 


rA-22rn 


99 


Expand  p,  u  and  u  in  the  power  series  of 
p  -  P  <°>  +  ap  <1>+— 


U  =  U 


(ti  =  Ui 


+  nu 


+  n<i) 


<1>+— . 
<1>+—  . 


(A-23a) 

(A-23b) 

(A-23c) 


We  confine  our  attention  to  the  limit  of 


0,  which 


corresponds  to  the  lowest  engenvalue  for  the  case  without 
swirl,  or,  in  the  terminology  of  Sozou  and  Swithenbank, 
to  a  slow  wave.  We  substitute  (A-23)  into  (A-21)  to  (A-22) 
and  retain  the  terms  up  to  0(n) .  This  yields  the  following 
expression  for 


i(1)  =  +  i\m\~  1  + 


(A- 24a) 


and,  a)  becomes  from  (A-23c) 


a)  ■  +  (|m|-  1  +  A_2lnil)n 


(A- 24b) 


where  +  sign  corresponds  to  m)0,  -  sign  to  m<0. 

For  the  special  case  of  m  =  1,  this  is  reduced  to 


1 

?  1 


which  is  identical  to  equation  (9)  of  Sozou  and  Swithenbank, 
which  appears  to  have  been  obtained  by  the  inspection  of  their 
numerically  computed  values  (see  their  Figure  8(a)).  In  their 
table  3,  they  also  list  the  results  of  computed  u;  in  the 
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following  Table  2 ,  we  compare  their  results  with  our 
analytical  formula  (A-24b)  side  by  side,  the  latter  being 
shown  in  parentheses.  The  column  left  blank  is  unavailable 
from  Sozon  and  Swithenbank.  On  the  whole,  our  expression 
appears  to  be  satisfactorily  close  to  the  numerical  results, as 
long  as  the  swirl  Mach  number  .defined  as  Ot^/Aq(0)  by  them  .remains 
less  than  0.5  or  so  .  Though  not  to  be  included  here,  if 
X  is  not  too  close  to  unity,  the  amplitude  of  u,  (A-23b) 
is  also  in  good  agreement  with  their  Figure  6. 
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"Family  tree"  of  streaming. 
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disturbance  spinning 
in  the  direction 
opposite  to  steady 
swirl 
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Figure  4(a)  Frequency  vs.  swirl  Mach  mimber 

X  ■  1.1 

solid  lines:  numerically  computed  values 
chain  lines:  analytical  results  from 
equation  (A- 8) 


streaming  onl) 


threshold  =  1.09 


streaming 
.  only 


steady 

swirl 


Figure  5(a)  Tangential  streaming  near  the  inner  cylinder,  ' 
showing  the  reversal  of  streaming  at  threshold 
swirl  for  the  first  radial  mode.  \  =  2.0 

K  =  0,  m=l;  temperature  at  outer  cylinder  wall, 
20°C,  Pr  =  0.71;  7  =  1.4  . 

solid  lines;  inner  wall  temperature  maintained 
at  its  steady  value. 

broken  lines;  thermally  insulated  inner  wall. 


threshold  M 


igure  6(b) 


without  unsteady  disturbance 


with  acoustic  streaming 


Figure  10.  Metamorphosis  of  total  d.c.  swirl  for  a 
Rankine  vortex. 
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FIGURE  CAPTIONS 


Figure  1.  Definition  sketch. 

Figure  2.  Lower  deck,  middle  deck  and  core. 

Figure  3.  "Family  tree"  of  streaming. 

Figure  4.  (a)  Frequency  vs.  swirl  Mach  number.  X=  1.1. 

Figure  4.  (b)  X  =  2.0. 

Figure  4.  (c)  X  =  5.0. 

Figure  4.  (d)  X  =  7.0. 

Figure  5.  (a)  Tangential  streaming  near  the  inner 

cylinder,  showing  the  reversal  of 
streaming  at  threshold  swirl  for  the 
first  radial  mode.X  =  2.0. 

Figure  5.  (b)  X  =  3.0. 

Figure  5.  (c)  X  =  5.0. 

Figure  5 .  (d)  X  =  7.0. 

Figure  6.  (a)  Threshold  swirl  Mach  number  for  the  first 

radial  mode  vs.  radius  ratio,  m  =  1.  Above 
the  threshold,  the  direction  of  the  tan¬ 
gential  streaming  is  opposite  to  that  of 
steady  swirl;  below  it,  the  direction  is 
the  same. 

Figure  6.  (b)  m  =  2. 

Figure  7.  (a)  Tangential  streaming  near  the  outer 

cylinder. X  =  2.0. 

Figure  7.  (b)  X  =  5.0. 

Figure  7.  (c)  X  =  7.0. 

Figure  8.  Metamorphosis  of  total  d.c.  swirl  for  a  free 
vortex. 


Figure  9.  Tangential  streaming  near  the  tube  periphery, 
showing  that  its  direction  is  always  the  same 
as  the  steady  swirl.  T/  (r.A  (r=r  ))  =  0.2; 

1  GX  O 

X  =  2;  k  *  0;  first  radial  mode;  P  =  0.71; 
y  “  1.4;  wall  temperature  is  alwayl  maintained 
at  its  steady  value,  20°C. 


